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We report on a fully nonequilibrium theory of the scanning 
tunneling microscopy (STM) through resonances induced by 
impurity atoms adsorbed on metal surfaces. The theory takes 
into account the effect of the tunneling current and finite bias 
on the system and is valid for arbitrary intra-adsorbate elec- 
tron correlation strength. It is thus applicable to the recent 
STM experiments on Kondo impurities. We discuss the finite 
temperature effects and the consequences of atomic scale res- 
olution of the STM for the spectral property of such systems. 
We find that the tip position affects the resonance lineshapes 
in two ways. As a function of distance from surface, the line- 
shapes vary due to the different extent of the adsorbate and 
metal wavefunctions into the vacuum. However, we do not 
expect large variations in lineshapes unless the tunneling into 
the tightly bound adsorbate states is considerable or when 
nonequilibrium effects are significant. As a function of lat- 
eral tip position, lineshapes should not change significantly 
on length scales of _R|| < 10 A under typical experimental 
conditions when the electrons tunnel into the perturbed bulk 
conduction states hybridized with the outer shell sp adsor- 
bate orbitals. Tunneling into surface states on (111) surfaces 
of noble metals should be important for observation of the res- 
onance at larger distances (> 10 A) and oscillatory variations 
in the lineshape should develop. This long range behavior has 
not been resolved in the recent experiments with Kondo impu- 
rities. The temperature dependence of the Kondo resonance 
cannot be deduced directly from the differential conductance 
as the thermal broadening of the tip Fermi surface produces 
qualitatively similar effects of comparable and larger magni- 
tude. Careful deconvolution is necessary to extract the tem- 
perate dependence of the Kondo resonance. The finite-bias 
current-induced nonequilibrium effects in tunneling through 
Kondo impurities should produce characteristic broadening 
of the resonance in the case of strong hybridization of the 
discrete state with the STM tip. 

PACS numbers:61.16.Ch,72.10.Fk 



I. INTRODUCTION 

A considerable body of experience and wisdom within 
the area of solid state tunneling phenomenon was built 
up throughout tunneling's "Golden Era of the Sixties" . 
It was during this period that many of the defining fun- 
damental ideas, basic theoretical strategies and method- 
ologies, and broad scope of new applications for tun- 
nel structures were first realized. A general introduc- 
tion to many of these achievements can be found in 



a number of comprehensive volumes^"^ and in the No- 
bel Lectures of Esaki, Giaver, and Josephson, who were 
awarded the 1973 Nobel Prize in Physics for "their [in- 
dependent] discoveries regarding tunneling phenomena 
in solids"**. It is against this background that the as- 
tounding achievements in contemporary tunneling stud- 
ies utilizing the single atom spatial resolution of the scan- 
ning tunneling microscope (STM) are most meaningfully 
considered^"**^. One phenomenon of key interest here 
which was first considered in the "Golden Era" is that 
of impurity/adsorbate-assisted elastic tunneling. Two 
bodies of work are particularly relevant to the present 
study. The first is the recognition by Appelbaum and 
coworkers of the possible role of the Kondo effect ^^^^'^ in 
determining certain current-voltage characteristics (e.g. 
"zero-bias anomalies") of metal-oxide-metal tunnel junc- 
tions containing localized paramagnetic impurity states 
near the metal-oxide interfaces'^. Second are the reso- 
nance tunneling studies involving valence electronic levels 
of single atoms adsorbed on metal surfaces, as probed in 
a field emission microscope configured for energy analysis 
(thus enabling electron spectroscopy) of the field emitted 
electrons'^"'^. Many years later, useful parallels between 
the theory of single atom resonance tunneling developed 
in the "Golden Era" and the theory of the STM, in the 
single-atom-tip limit, were unambiguously established'^. 
Further discussion of these issues from the past will be 
offered throughout the text, when appropriate. 

The basis for continuing interest and excitement in 
impurity/ adsorbate- assisted tunneling is that the trans- 
parency of tunnel junctions can be dramatically enhanced 
by the presence of states localized within the barrier 
when they are in resonance with the tunneling electrons. 
Tunneling through such states is, for example, the ori- 
gin of conductance fluctuations quantum dots exhibit in 
the Coulomb blockade'^. The tunneling probability in 
the presence of a "barrier" state is proportional to the 
spectral density produced by the hybridization of the lo- 
calized state with the conduction electrons and in many 
situations the current is given by the Breit-Wigner for- 
mula 
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Here F is the width at half maximum of the resonance 
produced by the hybridization with the conduction elec- 
trons in the right and left lead and eo is the energy of the 
local state. The value of F depends on the height and 
width of the barrier potential between the central region 
and the leads. 



Recently, enhancements in the zero bias conductance 
in quantum dots due to the Kondo effect have been ob- 
served^". It had been shown earher that {dl/dV), the 
zero-bias differential conductance, is proportional to the 
Fermi level density of states of the Kondo resonance on 
the quantum dot. Similarly, the Kondo resonance has 
been spectroscopically observed on single magnetic im- 
purities adsorbed on metal surfaces using the STM^^'^^. 
However in the case of the spectroscopic STM experi- 
ments, the resonance at the Fermi level appears to have 
an asymmetric shape and cannot be interpreted simply 
in terms of the local density of states of the impurity 
atom. Rather, the electron tunneling current being 
a coherent quantum effect - is a result of interference 
between competing tunneling channels, as will soon be 
detailed. Unlike in the quantum dot where the tunneling 
can take place with appreciable magnitude only through 
the quantum dot region, the apparent tunneling current 
from the STM tip to the surface can either go through 
the resonance localized on the impurity or directly into 
the conduction states of the surface. The distinction be- 
tween the conduction and local states will be discussed 
later. The notion that the tunneling conductance is pro- 
portional to the local density of states near the STM tip 
must then be modified. 

In addition to its most common use for observ- 
ing/determining atomic geometrical structure at sur- 
faces, the STM is used as a sensitive probe of surface 
electronic structure. Various theoretical approaches to 
the STM conductance employ the tunneling Hamiltonian 
introduced by Bardcen^"^ and Golden Rule type expres- 
sions in which under certain limiting conditions of prac- 
tical interest the STM conductance is indeed determined 
by the surface density of electronic states near the STM 
^jp24,25_ Xersoff and Hamann^^ developed a widely used 
model of the scanning tunneling microscope that includes 
the three dimensionality and spatial resolution of the tip. 

The generic problem of a discrete state interacting with 
a continuum of states arises in many different areas of 
physics and chemistry^^. In condensed matter physics a 
frequently occurring realization is the electronic state of 
an impurity atom immersed within a host lattice. In 
the case of magnetic impurities, the interaction gives rise 
to nontrivial phenomena such as the Kondo effect^ 

Within the context of atomic physics, Fano discussed 
related effects, as they might appear on observable ab- 
sorption lineshapes or resonant electron scattering cross 
sections which are due to the configuration interaction 
(CI) that couples a discrete two-electron excited atomic 
state with a continuum of ionization states^^''^". Whilst 
the "natural" lineshape of the resonance is Lorentzian, 
when studied by experiments in which an external probe 
interacts with the system, the resonance can appear to 
have an asymmetric lineshape. Such lineshapes are re- 
ferred to as Fano resonances. Fano found that an asym- 
metry in absorption lineshapes is due to interference 
between the excitation or decay into Cl-mixed discrete 
and continuum states which both couple to the external 



probe. If the coupling between the probe and contimmm 
is expressed in terms of an energy independent matrix el- 
ement tc and the interaction between the probe and the 
localized state (which has already been diluted by admix- 
ture into the continuum) by the matrix element ia, then 
the lineshape detected has the form 
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where q = ia/i^TrVtc) with V being the hybridization 
(or CI) matrix element between the local state and con- 
tinuum. The latter coupling results in the discrete state 
acquiring a width F = 27r PsV^ {ps is the density of con- 
tinuum states). 

In the present paper, we consider the problem of the 
discrete state embedded in a continuum using a probe 
such as the STM that has atomic scale spatial resolu- 
tion. This work has been motivated by the recent STM 
experiments involving single Kondo impurities^^'^^. The 
resonance observed in the conductance was interpreted 
by the authors in terms of the Fano interference. The 
fit of the resonance to the Fano formula^^ - generalized 
to the case where intra-atomic Coulomb interaction on 
the impurity is taken into account - was based on the as- 
sumption that the tunneling into an Anderson impurity 
can be extended to include the tunneling into the contin- 
uum in a straightforward way. Upon further considera- 
tion, it appears that this generalization of the Fano result 
to the case of STM conductance is not as straightforward 
as has been presumed. In the present paper we obtain a 
more complicated expression than the elementary Fano 
formula, one which accounts for the correct asymptotic 
behavior for large tip-impurity separation. In particular, 
when the dependence of the probe's distance from the 
local state is properly included, then observable conse- 
quences of the local state admixture with the conduction 
electrons show the correct asymptotic long range behav- 
ior in the largo tip-to-impurity-separation-limit as they 
obviously should. 

The difference between the Fano lineshape and the 
result obtained here is due to the different nature of 
the probe in the "multi-center" STM configuration (one 
"center" on the impurity /discrete state, the other on the 
STM tip/probe) compared with the "single-site" atomic 
physics processes. While Fano was concerned with light 
absorption or electron scattering whore the system un- 
der study was always at the probe's focus, the outcome 
of STM experiments must depend upon the variable spa- 
tial position of the probe (tip) with respect to the dis- 
crete state under investigation. Put another way, for the 
atomic physics applications considered by Fano, both the 
discrete state coupled to the continuum and also the ini- 
tially excited decaying or autoionizing state (the "probe" 
state in our STM language) are atomic states spatially 
localized at the same site by the same atomic central po- 
tential. In contrast, the "S" (=scanning!) in STM assures 
that the tip, hence initial excited state, can be indepen- 
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dcntly located with respect to the position of the "dis- 
crete state coupled to the continuum" and it is this extra 
degree of freedom that enriches the potential informa- 
tion content in STM lineshape analysis, but also requires 
a much more detailed theoretical treatment than merely 
fitting the atomic physics Fano lineshape, Eq. (1.2), to 
position-dependent STM spectra. This will be expanded 
upon in depth later. This realization demonstrates the 
importance and crucial need for considering the measure- 
ment process in quantum mechanical observations. 

While the structuring of this comprehensive paper is 
based on a logical development of the subject matter, it 
may be useful to present a roadmap of key points and re- 
sults to guide the casual as well as the dedicated reader. 
In section II, we introduce the model of the system and 
discuss our approximations of the tunneling matrix el- 
ements resulting in Eqs. (2.13), (2.14), and (2.17). In 
section III and the supporting appendices (A-C), we de- 
velop the general nonequilibrium theory of STM tunnel- 
ing current and conductance in the presence of an adsor- 
bate induced resonance. The more familiar equilibrium 
limit, asymptotically exact for large tip-surface separa- 
tion but pragmatically useful even for moderate ~ SA 
separations, is treated as a special case in IIIB. The re- 
lationship between the equilibrium tunneling resonance 
lineshape observable in STM experiments and the asym- 
metric Fano lineshape is established in III C (Eqs. 3.24 
and 3.25). The nonequilibrium contribution to the line- 
shape treated as a correction to the equilibrium limit is 
taken up in HID. The differential conductance is intro- 
duced in HIE. The crucial role of the substrate elec- 
tronic structure is explicitly considered in section IIIF 
where the substrate Green's function is evaluated for a 
jellium surface and in III G where surface states and real 
electronic structure effects are discussed qualitatively. 

In section IV, we illustrate the predictions and con- 
sequences of our theory on two models for the adsor- 
bate: nonmagnetic and Kondo. Using a jellium substrate 
model, we first discuss the common features of the two 
adsorbate models in terms of the non-interacting Ander- 
son model in IV A. Four different families of spectro- 
scopic lineshape variations are taken up: dependence on 
(1) the relative strength of the tip-to-adsorbate vs. tip- 
to-substrate tunneling and on the substrate electronic 
structure; (2) tip-surface separation; (3) lateral tip po- 
sition; (4) temperature due to Fermi level smearing. We 
explain why no variations in the resonance lineshape 
should be observed for lateral tip positions on length 
scales ~ (1 — 2) A - characteristic of the bulk kp - 
and only small lineshape variations should be expected 
for vertical tip variations (in experimentally accessible 
range). We leave the discussion of the tip and bias ef- 
fects on lineshapes to section IV B 3. 

Tunneling characteristics specific to Kondo systems 
are presented in IV B. We begin with conceptual issues 
in IV Bl. We then discuss the recent experiments on 
Co/Au{lll)^^ and Ce/Ag{lll)^'^ and relate our work to 
related theoretical papers^^'^^ ( IV B 2). In particular. 



we show that the temperature dependence of the Kondo 
resonance is not easily extractable from the temperature 
dependence of the differential conductance. We reiter- 
ate that the stability of the experimental lineshape with 
the tip position is to be expected. Variations in the line- 
shapes would, however, occur at larger distance due to 
tunneling into the surface state. The effect of adsorbate- 
tip hybridization and bias induced nonequilibrium on the 
current (conductance) vs. bias measurements in Kondo 
systems are dealt with in section IV B 3. Finally, an enu- 
meration of specific conclusions is offered in section V. 

II. MODEL AND APPROXIMATIONS 

Models of scanning tunneling microscopy are abun- 
dant in the literature of the last two decades and stan- 
dard texts exist We approach the problem as a 
nonequilibrium process and discuss the corrections to the 
TersofF-Hamann formulation^^. Our intent is to develop 
such theory under general and self-consistent assump- 
tions that accurately capture most of the qualitative as- 
pects involved and do so in a way that make extension 
to more realistic calculations formally straightforward. 
We focus on the tunneling through adsorbate resonances. 
Throughout this paper, we adopt the convention that 
the energies are measured with respect to the respective 
Fermi levels of the substrate and tip unless specified oth- 
erwise and seth = 1 . When the tip is biased we explicitly 
shift the tip energies. 

A. Model of the studied system 

We consider a system which consists of a clean metallic 
surface with a single impurity atom adsorbed on it. The 
STM will be used to study the system by means of tunnel- 
ing through a resonance produced by an electronic state 
of the impurity, such as the 5/ orbital of Ce/Ag{lll)^^ or 
3d orbital in Co/^?i(lll)^-'^. Unless otherwise noted, we 
place the origin of the coordinate system at a point on the 
surface of the metal directly below the adsorbate. This 
means that the position of the impurity is Rq = (0, 0, Zq). 
The system without the probe is described by the degen- 
erate Anderson Hamiltonian 

H,{Ro) = J2 ^oiRoKca + J2 UiRo)nana' + (2.1) 

a a>a' 

+ ^^kci^'^ka + J2 {^ka{Ro)ci^Ca + H.C.} . 
ka ka 

Here, eo is the energy of the impurity state tpa{r}, which 
we assume may be a multiplet of states described collec- 
tively by the quantum number a = (mo") . In the simplest 
case, a correspond to the spin a, but it may also include 
orbital degeneracy (m) in more complicated cases. In this 
paper, we discuss at most spin degenerate states with 
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a = a and N = 2 (degeneracy). We denote by cl the 
creation operator for this state. The Ck is the conduction 
band state energy - independent of a in the absence of 
magnetic field - with c^^ being the creation operator for 
the corresponding Bloch state with symmetry (spin) a 
common with the impurity state, and Vka is the matrix 
element for hybridization between the impurity and con- 
duction states. The second term in (2.1) corresponds to 
the intra-atomic Coidomb interaction between electrons 
in the impurity state tpa- 
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FIG. 1. schematic picture of the STM 

If the renormalized energy eo lies within the conduction 
band, the bound state broadens into a resonance which 
in the wide-band limit and with = 0, has a Lorentzian 

shape 



Po(w) 
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where 



27r (w - eo)2 + (r/2) 
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(2.3) 



with ps and V the assumed-energy-indcpcndcnt density 
of conduction states and hybridization matrix element 
from (2.1). In most of our later considerations however, 
we will retain an energy dependence in the model density 
of states, Psiijj). 



B. Interaction of the system with the STM tip 

When the STM tip is brought near the impurity, elec- 
trons can tunnel between the tip and the adsorbate state. 
This situation is expressed by adding an interaction term 
to the Hamiltonian 

Ho.t{Ru Ro) = {tap{Rt, Ro)cicp + H.c.} . (2.4) 

pa 



The tip states arc denoted by subscript p while the con- 
duction states of the metal by subscript k. The transfer 
matrix tap depends on the position of both the adsorbate 
Zo and the tip Rt = (-R|j, Zt). 

If the STM only coupled to the discrete state with 
transfer amplitude ta then the conductance would, in the 
wide-band limit {tap = ta, independent of p), be deter- 
mined by 
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and the conductance would thus be directly related to 
the impurity density of states. This is reminiscent of 
the defining characteristics from field emission resonance 
tunneling spectroscopy in which tunneling from the sub- 
strate to vacuum is disproportionately smaller than that 
from "good" adsorbates^^"^^'^'^. However, since in the 
STM geometry, tunneling directly between the tip and 
the metal surface can be comparable to (or in excess of) 
that between the tip and the impurity, the conductance 
exhibits a more complex behavior than that of a simple 
impurity local density of states. We take such processes 
into account through 

Hst{Rt) = J2 {tkp{Rt)clcp + H.c.} , (2.6) 

pk 

where the tunneling matrix element tkp depends on the 
position of the tip Rt- In the rest of this section, we will 
address the issue of the probe's effect on the system itself. 
This is particularly important when the tip is brought 
very close to the adsorbate so that the tunneling tap and 
tkp are comparable with Vka- In this case, the width F of 
the adsorbate resonance is no longer given by (2.3) but 
rather by 



T = 2nipsV^+Pttl) 



(2.7) 



The effect of the tip on the conduction states can also be 
important. However, the most important consequence 
of the tip perturbation comes when finite bias is applied 
across the tunnel junction. The system will be out of 
equilibrium and the problem of the system-probe must be 
approached sclf-consistently with the tip included in the 
system it is probing. The typical operational mode of the 
STM during imaging and spectroscopic measurements is 
such that the tip distance from the adsorbate is several 
atomic units larger than the adsorbate-surface separation 
and therefore |tap|,|ifep| ^ \Vka\- In this limiting case, 
the probe has no effect on the adsorbate-mctal complex 
other than as a source of hot tunneling electrons. This 
is a reasonable assumption as long as the tip separation 
is large enough to justify the approximation |fap|) l^fepl ^ 

\Vka\. 
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C. Model of the STM tip 

An important property of the tip is its spatial reso- 
lution, as discussed in great detail by niany^'-'^*^^'*'^^"^^. 
We will consider the tip to be well defined and terminated 
by a single atom through which the tunneling predomi- 
nantly takes place. This is the s-wave model of Tersoff 
and Hamann^'*, see Fig 1. The important features are 
the following: (a) the tip Hamiltonian is 

Ht^J^^P^^'P' (2-8) 
p 

where creates an electron in the state tpp{'r) with en- 
ergy e-p measured from the Fermi level of the tip epi] (b) 
when the tip is positioned near the surface, tunneling into 
and out of a state Vp can take place; (c) the states are 
filled up to the chemical potential ept controlled by the 
bias; (d) the tip states are characterized by a density of 
states which we denote by yOt(w) = ^ ^p)' and (e) 

the asymptotic form of the tip eigenstate tpp in the vac- 
uum region extending towards the metal surface is char- 
acterized by the atomic orbital of the apex atom. The 
wavefunction ^jjp can be found based on simple physical 
arguments without solving the complete problem. If 
is the work fmiction of the tip and Kt = •\/2m^(^7^-e^, 
then following Tersoff and Hamann, 

Vp(rlaiie'^'^?^^Pt^ (2.9) 

where R is the radius of curvature of the tip about its cen- 
ter which is located at the origin of this "tip-defining" 
coordinate system. While (2.9) represents an "s-wave 
tip" , more generally ipp would carry whatever symmetry 
was possessed by the relevant atomic orbital centered at 
the tip apex^°'^^'^'*. The wavefunction tail, controlled by 
K, depends on the bias and tip-surface separation. Both 
factors modify the height of the vacuum barrier, hence 
effective work function (f> determining k. These modifi- 
cation can be essentially included by renormalizing the 
wavehmction tails and the densities of states by position 
and energy dependent factors via the tunneling matrix el- 
ements. We discuss the tunneling matrix elements next. 

D. Approximations for the tunneUng matrix 
elements 

An important role in our formulation is played by 
the tunneling (hybridization) matrix elements Vka, tap, 
and tkp since they include the dependence on electronic 
structure and the tip and adsorbate position. The de- 
sired quantitative accuracy of the model for the tunnel- 
ing process is to a large degree determined by the ap- 
proximations made in the evaluation of these matrix el- 
ements. We do this first for a general adsorbate- 
metal system and then for a jellium model. It is rather 



straightforward to include band structure effects using a 
realistic electronic structure calculation of the substrate 
Green's function. 

We begin with the discussion of the matrix elements 
Vka and tkp that contain the metal wavefunctions tp^{r). 
They have the form 

Mki{Ri)= j d'rr^{r)vsi{r;R)il^'i{r), (2.10) 

where Vsi is the potential representing the mutual inter- 
action of the two systems. The wavefunction V'fc(^) is a 
Bloch state of the unperturbed metal and V[(^ is 
wavefunction (in the coordinate system of the metal) of 
the adsorbate state a in the case of Vka or the tip wave- 
function in the case of tkp- Either way, ^i';(r) can be 
written in terms of the wavefunction with the origin at 
the adsorbate (tip) as V'KO ~ ^i)^ where Ri is the 

position of the adsorbate [Rq) or the tip apex atom {Rt) 
measured from a reference point on the surface. The en- 
ergies Ek of the metal electrons are written in terms of the 
perpendicular and parallel components as Cfe = efc^ -I- ek\\ ■ 
We follow the convention that Ck^^ is measured from the 
bottom of the 2-D band and e^^ is measured with re- 
spect to eps- For example, we write for the jellium model 
Efez = kl/2ml - D and Cfen = fcjj/2m*, where (-£>) is the 
energy of the bottom of the band with respect to the 
Fermi level. The Bloch states can generally be written in 
the relevant region outside the metal as 

^nki.^ = e-'^"'^«„fe„ (P, ^)e''^^-'- (2.11) 

where n is the band index, u^^ (p, z) is a function weakly 
dependent on z outside the surface and periodic in p, the 
electron coordinate in the plane of the surface. This form 
is equally valid for the metal band gap surface states that 
seem ubiquitous to STM studies on (111) noble metal 
surfaces since both the z-propagating Bloch states and 
the localized surface states are eigenstates of the same 
Hamiltonian with different eigenvalues at a given k. At 
the surface, evanescent states into the bulk that appear 
upon analytic continuation of the band structure into 
the domain of complex k-vectors cannot be rejected on 
the basis of physical considerations as they were for the 
perfectly periodic interior of the solid. 

For jellium, u^^j:^^ (r) is constant and the the metal 
states are then simply plane waves along the surface with 
exponentially decaying amplitude into the vacuum. Here 

Kns = ^y2m*{(j)s - f-nkj = ^2m*(0s - e„fe + e„fe||) with 

m* the metal electron effective mass (number) and (j)s the 
height of the tunneling barrier for a Fermi level electron. 
For bias voltages much smaller than the work function, 
(pg is equal to the metal work function. We omit the 
band index n in the rest of the paper unless we explicitly 
discuss the electronic structure effects. We also define 
= y/2m*{(j)s — co), where the energy factor (0s — iv) 
represents the effective tunneling potential barrier for an 
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electron with energy lu. For small bias voltages, « 0, 
and we can replace A^^ , which depends weakly on energy 
in this range, by its Fermi level value (= A). 

We shift the integration variable in the integral (2.10) 
to an origin centered on the adatom (tip), which gives 



Jz>-Zi 



(2.12) 



ut{Ri) 



-Vsair,Ri)Mr)- 



Matrix elements of this type have been the focus of in- 
tense study in the context of charge transfer processes 
at surfaces.-" *^ The main contribution to the integral 
in (2.12) comes from the region just outside the surface 
where the z-dependence of the metal states is that of a 
decaying exponential. In the case of both the adatom 
and the tip, the integrand is assumed to be well localized 
to the z > —Zq region by virtue of the spatial properties 
of ipi. Consequently, the integral will be reasonably con- 
stant for all fc|| < l/n where ri, the radial length scale 
for the atom or tip, sets the range of fcy . Since, as will 
soon be demonstrated, other factors in the full problem 
will provide a much more severe fc|| ciitoff, it is sufficient 
to represent the integral by a constant. We then write 
the matrix element approximately Mki{Ri) — Motp^{Ri), 
where Mq is the overlap integral defined in (2.12) that 
contains the dependence of the matrix element on the 
symmetry of the atomic orbitals near the tip and on se- 
lection rules. 

Since the distance of the adsorbate from the surface is 
small, it will be sufficient for the purpose of Vka to write 
the Bloch state in the form Vfe(^o) = c-^o/a^^^q) with 
the decay constant independent of the k- vector. With the 
choice of our coordinate system Rq = {0,Zq), we write 
Vka in the form of a separable product of a k-independent 
function of adsorbate position multiplied by a function of 
k 



VkaiRo) ^ Va{Zo)rM 



(2.13) 



where Va{Zo) = Voe~^°^^. We cannot make this last 
approximation in the tip-to-surface tunneling matrix el- 
ement tkp, since the interplay between k and Zf depen- 
dence has an important role in the tunneling process. 
However, from the conceptual point of view, we find it 
convenient to isolate a fc-independent dependence on the 
tip position in tc{Zt) = foe~^*/^ and write the matrix 
element in the form 



tkp{Rt)=tc{Zt)e^*/^rk{Rt 



(2.14) 



We note that the s-wave tip-to-surfacc tunneling matrix 
element tkp has been shown by Tersoff and Hamann^^ un- 
der quite general assumptions to assume the form (2.14) 
with the tip wavefunction given by Eq. (2.9), which indi- 
cates that our simple qualitative arguments seem to be 
supported by more detailed analysis. 



The tip-to-impurity matrix element 

tap{Rti Rq) = 



(2.15) 



/ 



(fr V;-(r - Rt)vat{r; Rt, Ro)'^a{r- Rq). 



depends on the position of both the adsorbate (Zg) and 
the tip {Rt)- Since the tunneling from the tip takes place 
predominantly through the apex atom, the wavefunction 
ipp in the last expression, a generalization of (2.9), can 
be written (with the tip at origin) as 



V'p-(r) = e 



-V's(r) 



(2.16) 



where Vs is the angular part of the orbital localized on 

the apex atom whose center is located at Rt . The tunnel- 
ing matrix element tap will thus reflect the symmetry and 
spatial dependence of the states localized on the adsor- 
bate and in the tip apex. Within the s-wave tip model, 
t/'s is .just an innocuous constant. Since '0p and V'a ap- 
pearing in Eq. (2.15) are both atomic-like functions in 
the relevant region of overlap centered respectively on 
the tip and on the adatom, in a broad sense tap is similar 
to a common two-center hybridization/hopping integral 
deflning the binding in a diatomic molecule^'^'^'*. Typi- 
cally the magnitudes of these integrals are exponentially 
decreasing functions of their separation (possibly multi- 
plied by a mildly oscillatory function accounting for the 
nodal structure of the atomic functions). Based on this 
analogy from quantum chemistry, tap given by Eq. (2.15) 
should take the form 



tap{Rt,2o) « tae' 



-\Rt-Zoi^\/a 



= ta{Rt,Zo) (2.17) 



Here, a ^ w (^Kt + -^2m(^s — eo)^ is an effective decay 
constant evaluated f or states at the Fermi level of the 
tip, |_R( — Zoi^l = \ R\ + {Zt — Zt^Y is the tip-to-atom 



separation, as depicted in Fig. 1, where R\\ is the parallel 

component of Rt- The decay constant Kt depends on the 
energy of the tip state (cp), but this dependence is very 
weak for small biases considered here and we neglect it. 
In this case, the matrix element is well approximated by 
the p independent form ta{Rt, Zq)- The matrix element 
ta may be taken real. 



III. THE NONEQUILIBRIUM THEORY OF THE 
TUNNELING CURRENT AND DIFFERENTIAL 
CONDUCTANCE 

The tunneling between the adsorbate-metal complex 
and a biased tip is a nonequilibrium process. Although 
we frequently make the assumption in this paper that 
the tip-system interaction is weak enough so that local 
equilibrium is maintained to a good approximation, the 
assumption is less valid when the tip is near the surface. 



6 



For this reason, we develop our theory of the tunnel- 
ing process within the Keldysh-Kadanoff'*^'''^ framework 
for the nonequilibrium Green's functions and discuss the 
nonequilibrium corrections. 



A. General expression for the tunneling current in 
terms of nonequilibrium Green's functions 

We define the tunneling current as the flow of electrons 
through a closed surface around the tip. It is expressed 
in terms of the continuity equation as 



dt 



(3.1) 



where 



is the number operator for the tip 



electrons, and the brackets signify the ensemble average, 
which in the local equilibrium case is the thermal average 
over the tip states. The time derivative is found from the 
Schrodinger equation of the total Hamiltonian of the tip- 
substratc-adsorbate system Hioi = Hs + Hi + Hai + Hsi- 
Since the number operator commutes with Hs and Ht, 
the only contribution comes from the interaction terms 
Hat and Hgt and the current is 



7 = ^ Im {^tkp{ 

kp 



ap 



where we omitted the arguments in tkp and tap for con- 
venience and used the relation tap = t*a and tkp = t*/^. 
We write the arguments explicitly only when we wish 
to emphasize their dependence. We define the time 
loop Green's functions Gpa{t,t') = —i{TcCp{t)cl{t')) and 

Gpk{t, t') — —i{TcCp{t)c\{t')), where Tc orders the times 
along a contour C in the complex time plane. The con- 
tour can be taken to be the Kadanoff-Baym contour^^, 
the Keldysh contour^^, or a more general choice. The 
discussion of nonequilibrium Green's function is available 
in standard books and review articles'*^"^* and we refer 
the reader to these references for further details. Our 
notation follows closely that of the more detailed discus- 
sion in reference^®. In the present paper, we are mainly 
interested in steady state tunneling current (time inde- 
pendent) and therefore we work with the Fourier trans- 
formed quantities in frequency rather than time space. 
The current can be written as 



O poo 

7 = - Im / du^iY, tkpG<,{co) + J2 tapG<a{^)}, 
" kp ap 

(3.3) 

where G^j.{lj) and Gp„(a;) are the Fourier transforms of 

GtS^n = (4(i')cp(t)) and G<a{t,t') = (ct(t')cp(t)) - 
the analytic pieces on the real time axis of the Green's 
functions introduced above. Equivalently, the current 



may be calculated from + ^^)- "^asy to see 

that this approach also leads to the equation (3.2). 

The problem of finding the current thus reduces to 
finding the "lesser" Green's functions G^a and G'^j,. This 
is done using the equation of motion method for the time 
ordered Green's functions in Appendix B and the rules 
for analytic continuation described in Appendix G. In 
order to see the interference between the two scattering 
channels giving rise to the Fano lineshape, it is also useful 
to express the current using (B2) and (B6) in Eq. (3.3) 
as 



7= — Im 
h 



/oo 
d. 
p -oo 

{G^^^ \tpa\^ Ga-\- ^ taptpk Gka + 

a ka 

H" ^ ^ ^kpipa Gak H~ ^ ^ ^pk' ^kp ^/c'/c]} • 



(3.4) 



ka 



kk' 



The first term corresponds to the tunneling into the ad- 

sorbate state hybridized with the metal electrons. The 
fourth term is the contribution to the current from the 
direct tunneling into the conduction states perturbed by 
the presence of the discrete adsorbate state. The second 
and third terms give the interference between the two 
channels. 

We substitute the solutions (B4) and (B9) for Gpk and 
Gpa from Appendix B into Eq. (3.3) and write the tun- 



(3.2) neling current in the form 



7= ^Im 

h 



(3.5) 



X ^ {G^^,{^)T<^{co) + G<^,{^)tX{co)} , 
pp' 

where the retarded function T^,{Rt,Ro, to) plays the role 
of the T-matrix for scattering of the tip electrons from 
the adsorbate-metal complex. It is defined by 

^pp' ~ 5Z ^pk^kikp' + ("^-6) 

fc 

+ ipaGa {tap' + VakG^t^p'), 



with 



and 



Vka 



ipa — ^pa H~ ^ ^ ^pk G}^ Vfco, 



(3.7) 



(3.8) 



the hybridization and tunneling matrix elements for the 
adsorbate modified by the tip-substrate interaction. The 
matrix Tpp/ incorporates the properties of the tip as well 
as the adsorbate into the expression for current. We dis- 
cuss its physical meaning more in the next section. 
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B. Equilibrium limit of the tunneling current at 
large tip-surface separation 

We define the equilibrium tunncUng current as the 
large tip-surface separation limit of (3.5) when the tip, 
adsorbate, and substrate are all in local equilibrium. 
This is equivalent to keeping only the lowest order terms 
in tkp and tpa- In our formalism, this is achieved by 
replacing Gpp' G°Spp' in Eq. (3.5), Vka Vka in 
Tpp' , and by using the fluctuation-dissipation relation 
Gf{u!) = fi{uj)pi{uj). The subscript i stands for tip (t), 
adsorbate (a), and metal (s), respectively. The adsorbate 
is in equilibrium with the metal, i.e. /a('A') = fs{<^)- The 
matrix Tpp' is expressed entirely in terms of the Green's 
functions of the system and the tunneling matrix ele- 
ments tpk, tpa- Since these matrix elements reflect the 
symmetry of the apex atom wavefunction, but are only 
weakly dependent on p on the energy scale of the res- 
onance width, the matrix Tpp> will also have this prop- 
erty. We therefore make an additional assumption that 
EpGOTpp~(EpG°)T*, where 

Tt = tpkG^tkp + ipaGaiap (3-9) 
k a 

is only a function of the atomic tip orbital independent 
of p. We define a tip-specific quantity observable by the 
STM, which is related to the local density of states 

p,at{Rt,Ro;to) = --ImT,«(i?t,i?o;'^), (3.io) 

TT 

and write the equilibrium current Ieq{Rt, Ro, V) as 

Ieq = -^J dui[ft{ui') - fs{ui)]pt{ui')psat{ui), (3.11) 

where pt = J2p^i^ ~ ^p) ^hc density of tip states and 
u)' = w — eV with V being the bias voltage. This equa- 
tions is easily related to traditional formulations given in 
terms of an integral product of an electron "supply func- 
tion" multiplied by a tunneling or transmission probabil- 
ity^~^'^^'^^ when it is realized that psat, as defined here, 
already contains within it factors (ex |t|^) representing 
the role of the tunneling probability. 

This expression has a form similar to the standard tun- 
neling theories which express the current as a product of 
the local densities of states of the two systems evaluated 
at a common point and a difference in the correspond- 
ing Fermi functions. Kawasaka et al.^°'^^ who studied the 
STM current through a Kondo resonance used as a start- 
ing point of their considerations the Tersoff and Hamann 
expression^^ 

IeqCxJ(Lj [ftiij') - fsiij)] pt{uj')psa{RuOj) (3.12) 

according to which the current at zero temperature is re- 
lated to the LDOS of the adsorbate plus metal electrons 



Psa{Rt,^^) at the position of the tip, where the local den- 
sity of states is 

PsaiRuco) = --Im {Rt \G''{co)\Rt). (3.13) 

TT 

The LDOS is expressed in terms of unperturbed metal 

and adsorbate states by inserting \k){k\ + |a)(a|(w 1) 
on both sides of G in (3.13). This is strictly valid only for 
orthogonal orbitals, {a\k) = 0. The four resulting terms, 
proportional to Ga, Gka, Gak, Gkk', reflect the fact that 
the LDOS includes both the adsorbate and metal states 
perturbed by their mutual interaction. Inserting this ex- 
pansion into Eq. (3.13) gives for psa{Rt,i^) 

p,a = ~Im {J2 + ^ ^j: + (3.14) 

a ka 
ka kk' 

where the wavefunctions are evaluated at Rt- The four 
terms correspond to the terms in the square brackets 
in (3.4). We can also rewrite psa using the expres- 
sions for Gak and Gkk' in Appendix B. We deflne 
V'a = V'a + Efc i'kGlVka and write 

Psa = -^Im \^k? G'^ + J2 l^aP Gt}. (3.15) 

fe a 

Comparison of (3.15, 3.12) with ( 3.9, 3.11) shows the 

difference between our equilibrium limit and the trans- 
fer Hamiltonian method^'^^. The tunneling current and 
differential conductance in the equilibrium limit provides 
information about /7.,o/ a local density of states modi- 
fied by the tunneling matrix elements rather than the 
LDOS. In the case when the tunneling takes place into 
distinct orbitals with different symmetry, psat can be 
rather different from psa and the statement that the STM 
is a measure of local density of states must be understood 
in this context. 



C. The equilibrium tunneling current and the Fano 
lineshape 

In this section, we evaluate psat in Eq. (3.11) and write 
the equilibrium current using the approximations (2.13), 
(2.14), and (2.17) for the tunneling matrix elements. We 
note that the approximations do not require any speci- 
fication of the substrate electronic structure. The final 
form of the current allows the discussion of the tunneling 
resonances in terms of the well established Fano line- 
shapes. 

We introduce new quantities in terms of which the cur- 
rent is expressed. First, we define the "bulk" density-of- 
states (DOS) for the substrate and the tip as Ps{oj) = 
Y,k S{u; - ek) and pt(w) = Y^p S{uj - Cp). The impurity 
width without the STM tip is defined as 
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r„,(i?o,a;) = 27rp,(a;)y„2(iJo). 



(3.16) 



The adsorbate perturbation on the local density of con- 
duction states at some lateral position between the tip 
and the adsorbate is discussed in terms of the unper- 
turbed substrate Green's function 



G+(f,r";u;) = ^ 



LO - ek + irj 



(3.17) 



We define two dimensionless quantities related to the real 
and imaginary parts of the Green's function 



A{R,u)) = e 



zf^ Re G+{R, Q;uj) 



and 



^{R,uj) = -e 



irpsiuj) 



(3.18) 



(3.19) 



These two functions carry the information about both 
the spatial extent of the metal electron perturbation at 
arbitrary R in the surface region due to a localized per- 
turbation at _R|| = and also the spatial resolution of the 
tip, as we will sec later. We have included the exponen- 
tial factor e^/^ in the definition (3.18), (3.19), and (3.20) 
because we explicitly take the ^-independent part of the 
exponential dependence on position to be part of the tun- 
neling matrix elements Va{Zo), ta{Rt, Zq), and tc{Zt). 
We postpone further discussion of GQ{r,r' ,io) to the 
subsection IIIF. 

Finally, wc define a dimensionless quantity as the nor- 
malized density of the substrate states at a position R 
above the metal surface 



u{R,(jj) = e 



2Z/X 



„2Z/A 



Im G(j"(i?, R;uj) 



(3.20) 



The tunneling current Iq into a clean metal is given by 

the first term in Eq. (3.9). The current Iq for small bias 
{V/(j)s <C 1) can be written with the above definitions in 
a familiar form 



Io{Rt,V) 



2e r 



duj X 



(3.21) 



Xpt{uj') [ftiu') - f,{Lj)] p,{ij)TTtl{RtMRt;Lj). 

Here, /s(w) and /t(w) are the substrate and STM tip 
Fermi functions, respectively, and ui' = w — eV. The tip 
and substrate are assumed to have common chemical po- 
tential eps = ^Ft = at zero bias eV = and wc adopt 
the convention of measuring the energies in the substrate- 
adsorbate complex and in the tip from their respective 
Fermi levels at finite bias. The bias V is measured with 
respect to eps and is defined as positive when the chem- 
ical potential of the tip ept is raised. The functions ps 



and pt are the substrate and tip densities of states, re- 
spectively. 

It follows from Eq. (3.21) that the tunneling current /q 

is independent of temperature if pt, Ps, and v{Rt,Lij) are 
independent of energy in the relevant energy range. If ps 
shows structure on the scale of the temperature T while 
pt is constant, the current will depend on the temperature 
of the tip only, and vice versa. These statements are not 
limited to the case of clean metal surfaces, but also hold 
for the substrate with an impurity. The same is true for 
Q, the differential conductance. 

The equilibrium current in the presence of the adsor- 
bate is written by expressing psat with the notation and 
approximations that lead to Eq. (3.21). We define a mod- 
ified matrix element ta{Rt, Ro',u)) for tunneling from tip 
to the adsorbate state as 



ia = ta + TTtc Ps AVa- 



(3.22) 



The second term represents a coherent process of tip- 
to-surface tunneling, through-surface-propagation, and 
surface-to-adsorbate hopping. This is completely isomor- 
phic with Fano's coupling of an excited state (here the 
tip state) with the originally discrete state "modified by 
admixture of states of the continuum". The reader is 
enthusiastically directed to the original Fano paper for 
further enlightenment on this point. 

We introduce the Fano^^ parameter q{Rt, Ro;co) as 



t„ 



T^tcVaps 



(3.23) 



We will see later that this definition of q makes the 
expression for differential conductance formally equiv- 
alent with the Fano formula in certain limits. It is 
rather straightforward now to evaluate psat using (3.9) 
and (3.10) in (3.11). After rearranging the terms, we 
write the current Ieq{Rt, Ro,uj) in the presence of the 
adsorbate resonance as 

^ ^ 2e 1"°° 

Ieq{Rt,Ro, ^) = ^ j Pt^^') X (3-24) 

X [ft{oj') - fs{oj)]ps{oj) TTtliRt) Y{RuRo,u;), 

with 

^ = ^ + E ^ {(7' - q') Im + 2qj Re } . (3.25) 

a 

In our approximation, the localized nature of the tip 
and the adsorbate enters through the position depen- 
dence of taiRtjZo) and the substrate Green's function 
GQ{Rt,0;uj). The matrix element gives an exponen- 
tially decreasing amplitude with increasing tip-adsorbate 
distance and the substrate Green's function gives de- 
creasing amplitude due to the phase difference between 
electrons entering (or leaving) the surface at the adsor- 
bate site and leaving (or entering) at (.R|| ,z = 0) and also 
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due to the exponential decay of the tip wavcfunction with 
increasing fcy). We note that, in the wide band hmit for 
the substrate and with the tip near the surface above the 
adsorbate, ia ^ta, since in this hmit Re G+ and thus A 
vanish. 



(a) 



Equilibrium 



(b) 
0.16 eV 




(c) 



(d) 
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E (eV) E (eV) 

FIG. 2. The spectral function Pai^j) and the occupied den- 
sity of states of a model Kondo system as a function of bias 
voltage at small metal-tip separation given by Tat = O.lFas. 
(a) equilibrium, (b) - (d) finite bias. Bold line = electron 
population, solid thin = spectral density, dotted = equilib- 
rium spectral density. 



D. Nonequilibrium effects at stronger tip-surface 
coupling 

We now generalize Eq. (3.24) for the equilibrium tun- 
neling current - obtained in the lowest order in tap and 
tkp ~ by including nonequilibrium effects. The general 
problem of tunneling for arbitrary relative strength be- 
tween the tunneling amphtudes tap and tkp and the hy- 
bridization matrix Vak and for finite bias is formulated 
in Eq. (3.5), but the expression is quite complicated to 
evaluate in practice. In a typical STM experiment, the 
tunneling matrix elements tap and tkp are much smaller 
than than Vak ■ We can expect the nonequilibrium effects 
to be important when, at small separations, the magni- 
tude of the two tunneling matrix elements is not a neg- 
ligible fraction of |T4fe|. However, we can always safely 
assume that |tafe|, \tap\ are smaller than \ Vak\ in the STM 
experiments under all realistic conditions. 

Therefore we make additional simplifications which are 
justified by these relations. First of all. we replace Vak 
by Vak inside Eq. (3.6) and (A5). We neglect the modi- 
fications to the tip and substrate wavefunctions, i.e. re- 
place Gkk' by Skk'G^ and Gpp' by 5pp/Gp. We also ne- 
glect any deviations from thermal electronic distribution 
in the substrate and tip, i.e. we assume the validity of the 



fiuctuation-dissipation theorem for the tip and substrate 
Green's functions. On the other hand, when the tip- 
adsorbate coupling is not negligible with respect to the 
adsorbate-metal hybridization, the current into the reso- 
nance can be large enough to produce significant nonequi- 
librium electronic population on the adsorbate since the 
time scales for electron dissipation from the resonance 
into the metal and tip, respectively, are comparable. In 
this case, the fluctuation-dissipation theorem = fsPa 
is no longer valid for the adsorbate Green's function and 
we must use the full nonequilibrium G^{uj) instead of 
fs{^)Pa{^) in Eq. 3.24. 

Under these assumptions, we flnd it convenient to write 
the total current with the nonoquilibrimn effects as Itot = 
leg + SInon, where leq is formally given by (3.24) and 

^■^non is 

rj poo 

SInon = -- E / ^^PtpltlVo X (3.26) 

a J-oo 

X {fslmG^ + nG<){q^+^% 

where all adsorbate and substrate densities and Green's 
functions are evaluated at energy uj and pt at oj' = 
oj — eV. We omitted the spatial arguments for sim- 
plicity. The bias dependence enters through the self- 
consistent solution of the adsorbate spectral density 
Pa{u}) = — ilmG^(w) and the "lesser" Green's function 
G^(cj). In the case of nonintcracting system, (U = 0), 
the spectral density does not depend on the bias and the 
only nonequilibrium (finite bias) effect is given by the dif- 
ference between the equilibrium G'^ ^^{uj) = fa{ijj)pa{'^) 
and the nonequilibrium density of occupied states G^ (w) , 
as featured in 8Inon- 

On the other hand, the spectral density Pa(w) of Kondo 
systems depends on the bias. This means that leq also 
contains nonequilibrium effects and is different from the 
equilibrium current despite the subscript "eq" and its 
identical form. The effect of bias on the spectral func- 
tion depends on the tip hybridization with the discrete 
impurity level and is similar to that of temperature for 
eV < Tk where it broadens the Kondo resonance. At 
larger biases the broadening increases further and a sec- 
ond peak may develop at the Fermi level of the tip, 
depending on the strength of the adsorbate-to-tip hy- 
bridization Tat = '^'^Ptta compared to Tas = '^i^PsVa for 
the relevant impurity orbital. ^^'^^ In Fig. 2, we show for 
different bias voltages the spectral function and electron 
occupation of the resonance for a model Kondo system 
with Tat equal to ~ 10% of Tas and under an additional 
assumption that \tkp\ -C \tap\- The model will be dis- 
cussed in more detail in section IV B. 



E. Differential conductance in the limit of large 
tip-surface separation 

The differential conductance is obtained directly 
from (3.24) by differentiating it with respect to the bias. 
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i.e. Q = dl/dV. Wc do this here under the assump- 
tion that the bias vohagc is varied across a sufRcicntly 
narrow range so that the density of tip states may be 
taken constant. Under these assumptions the differential 
conductance Qeq is 

X p,{tj) TTtliRt) Y{Rt,Ro;w) 
and for the clean metal 

eo(4, V) = ^p cLoptiuj') (-^^) X (3.28) 
X ps{ij) irtliRt) u{Rt,ij). 

These expressions neglect any changes to the tunneling 
barrier from the finite bias voltage. When these approx- 
imation are not justified, the conductance must be ob- 
tained by differentiating the expression for current (3.24) 
and (3.21). This is always the case for 1^^^ of the previous 
section when nonequilibrium effects are important. 

It is known from the Anderson "compensation theo- 
rem" that in the wide band limit {ps{^) is constant and 
unbounded), the presence of an impurity does not affect 
the conduction electron density of states at all. However, 
the density of conduction states is affected locally even in 
this hmit, as can be seen by setting i/'a = in Eq. (3.14). 
The perturbation of the conduction electrons is probed 
by the STM directly if ta = 0, i.e. when tunneling into 
the local state is absent due to cither large tip-adsorbatc 
separation or because of symmetry. Under these condi- 
tions if the STM conductance shows the signature of the 
local resonance, it is a result of the perturbation of the 
LDOS of conduction electrons. 

An important feature of the final result (3.24) 
and (3.27) is that the role of the impurity resonance on 
the tunneling conductance is contained in the Green's 
function Ga of the local adsorbate state. It is then pos- 
sible to separate the problem into two steps. First, the 
solutions for the adsorbate and substrate Green's func- 
tions arc found for a given system; then the tunneling 
conductance is calculated using the solution in the ex- 
pression (3.24) or (3.27). Since the equations (3.24) 
and (3.27) have been obtained under very general as- 
sumptions, they can be used as a starting point in the 
study of a variety of tunneling problems with appropriate 
approximations for the Green's functions and the tunnel- 
ing matrix elements. We demonstrate this in the next 
section, where we first study noninteracting and then 
Kondo systems. If the approximations to the tunneling 
matrix elements employed here are too crude, the more 
general expression (3.5) or (3.11) must be used. 

We find that the lineshape [Q vs. V) of the adsorbate 
resonance depends sensitively on (1) the relative strength 
of the tunneling matrix elements tap and tfcp, (2) the 
perturbation of the conduction electron states, (3) the 
lineshape of the local resonance, i.e. the spectral function 



— ifm Ga and Re Ga, and on (4) temperature in the 
case when T ^ F, the width of the resonance. Since the 
observed perturbation of metal states and the tunneling 
matrix elements depend on the tip position, so will the 
lineshape. 

F. The substrate Green's function and 
perturbation of the conduction electrons: jeUium 
surface 

There arc two ways in which the adsorbate state affects 
the tunneling conductance: (A) direct tunneling into the 
discrete state; (B) perturbation of the conduction elec- 
tron states by the discrete state which consequently con- 
tributes to the tip-to-continuum tunneling current. Both 
contributions drop-off with increasing tip-adsorbate sep- 
aration. The direct tunneling into the resonance is con- 
trolled by tap{Rt, Zq) which is a function of the overlap 
between the tip and adsorbate wavefunction and thus de- 
cays exponentially with the distance. The perturbation 
of the continuum also vanishes at large distances from 
the adsorbate. However, its spatial extent shows a more 
complicated behavior and depends on the details of the 
electronic structure of the substrate in resonance with 
the broadened discrete state. It is anticipated that this 
contribution will show a significantly longer range infiu- 
ence than the decaying exponential, much in the spirit of 
Fricdcl oscillations. 

The position dependence of the perturbation enters 
through the Green's function GQ{Rt,0; uj). We note that 
the imaginary part 7 appears explicitly in the expression 
for conductance, (3.27), while the real part A enters the 
definition of ia- The STM-observable effects of the spa- 
tially dependent perturbation of the conduction electrons 
caused by the local impurity state are thus also controlled 
by (3.17). In our formulation, a non-trivial Gq is 
intimately related to the dependence of the tunneling el- 
ement tkp{Rt) on the lateral tip position ^||. With this 
in mind we focus on the specific problem of tunneling as 
a function of the tip-adsorbate separation. 

We consider a simple approximation for Gq" based on 
the assumption that in the relevant surface region the 
surface corrugations are smoothed out (jellium model) 
and both the Bloch and/or surface state '^/j (2.11) is given 
by 

ip^{f)o^ ^''r p. (3.29) 

The states with the smallest Kg have the longest tail into 
the vacuum region and thus will be the most important 
ones in the tunneling process. These are the states with 
the smallest Ck^^- It is then reasonable to represent Ks 
in terms of the Taylor expansion around the minimum 
of e^ii with efe equal to the bias. In most cases, it is 
reasonable to replace eu by its Fermi level value. We 
expand around its minimum as ek« ~ kfJ2m*, and 
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write K, 



A- 



Afcjj/2 plus higher order terms which 
we neglect. For states far from ep^ and at small tip- 
surface separation, the expansion should be made around 
a different value of A. For the purpose of this paper, it 
is sufficient to consider the Fermi level value A. We then 
write 



Mr) 



(3.30) 



Zt ^ oo and "f{Rt,Lu) is independent of the lateral tip 
position. Since at the same time A ^ and ta 0, 
the STM has no spatial resolution in this limit. As the 
tip moves closer to the surface the spatial resolution in- 
creases. In the limit Zt = 0, the integral in (3.32) can 
be evaluated and j{Rt,u}) = jo{kuiR\\) where jo is the 
spherical Bessel function of zeroth order. 



As we will show later, the second exponential e~^^*'^i|/^ is 
a measure of the tunneling current carrying k\\, the prop- 
erty that gives the STM tip its spatial resolution, and the 

third exponential, e*''!! '^H , controls the dependence of the 
tunneling current on the lateral tip position. 

With this approximation for Bloch states in the surface 
region the substrate Green's function (3.17) is 



G+iRt,0;io) 



k 



-AZtfe^/2gife||.fi|| 



\M0)\' 



u> - ek+iri 



(3.31) 



For the bulk band state propagation, it is easy to show 
using Eq. (3.19) that 

j{Rt,cj) = /'(ixJo(A^i?||Vl-a;2)e-^^*'='(i-"')/2 
Jo 

(3.32) 



and 



V 

7rps(w) JO 



2D 



deps{e) 



(3.33) 



where Jq is the zeroth order Bcsscl function and k^^ is 
the wavevcctor of the substrate state of energy oj. The 
normalized density of (STM-accessible) conduction states 
viZt) a distance Zt from the surface is 



p{Zt,uj) = / dxe 
Jo 



-\Ztkl{l-x^) 



(3.34) 



In calculating A, 7, and v we assumed jellium-likc dis- 
persion relation uj — k'^/2m* and use parabolic density 
of states Psi^) = 1 — w^/i?^. The incompatibility of the 
density of states with the dispersion relation is not im- 
portant for the purpose of demonstrating the important 
band struc;ture effects at this level of simplification. 

Although the expressions (3.32)-(3.34) are valid for a 
very simple model of the surface, we believe they contain 
the most important features of more realistic bulk elec- 
tronic structures. We now discuss these features begin- 
ning with 7, Eq. (3.32). At large Zt, the dominant contri- 
bution to the integral in 7 comes from small values of the 
argument y in Jo{y)- In this case, we use the mean value 
theorem to write (3.32) as j{Rt,ui) = Jo(fci?||)e-^^'^'/2 
where k = aka, with a G (0,1). Clearly, a — > as 
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FIG. 3. (A) and (B) show the typical spatial dependence 

of the tunneling matrix element ta (bold) and of the sub- 
strate Green's function for the jellium model evaluated at 
= 1.2 (left panels) and k^ = 0.6 (right panels) for 
different Zt. We parameterize ta oc e"^!!''" with a = 0.75 A. 
Also shown is the Bessel function Jo(fca;-R||) for the same kui, 
i.e. 7 of a surface state with k^. The lower panels (C) and 
(D) show A, the upper panels (A) and (B) show 7, Jo, and 
ta. The comparison for the two wavevectors assumes identical 
barrier (damping constant A). 

The resonance lineshape depends on the ratio of the 
direct tip-adsorbate tunneling amplitude {ta) to the am- 
plitude for tunneling into the perturbed metal states (tc) 
and on the interference of conduction electrons scatter- 
ing from the impurity. The latter contribution is repre- 
sented here by the substrate Green's function G+ . There- 
fore, the presence of an impurity on the surface can be 
sensed spcctroscopically even if the direct tunneling into 
the resonance is negligible as is the case, for instance, of 
Ce/Ag{lU)^^. As our estimates for ta{R\\) and 7(i?t) 
indicate, the direct timncling matrix element ta{R\\) falls 
off much more rapidly with than does 7. This is due to 
the limited spatial extent of the tightly bound impurity 
orbital. Therefore, the relative importance of tunneling 
into the perturbed continuum is likely to increase with 
the lateral tip-adsorbate separation. 

Using the simple model for G+, we show the typical 
length scales in Fig. 3(A). We plot ta (bold solid line) 
parameterized as ta oc e^^^i^" and normalized to one 
for = together with 7 evaluated at three different 
positions Zt above the surface and with k^ = 1.2 A~^. 
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The exponential fall-off for conduction states at the Fermi 
level with work functions in the range 4 — 5 eV would be 
characterized by A = 0.9 A. The tightly bound discrete 
state will have larger decay constant and we parameter- 
ize it by a = 0.75 A. Clearly, 7 decays much slower than 
ta at small tip-surface separations, but the difference in 
fall-off becomes smaller with increasing Zt . We also see 
by comparison with the Bessel function Jo{kpR\\) (light 
dotted) that the spatial frequency decreases with increas- 
ing Zt and the oscillations eventually disappear entirely. 
This is due to the fact that smaller fc-vectors have larger 
weight in the tunneling at greater Zt (see integral (3.32)). 
The real part A, shown in panel (C) for the same k-vector, 
has a similar behavior. However, being a Hilbert trans- 
form of the imaginary part, the nodes in A appear at 
the positions of local extrema of 7 and vice versa. As 
we will see later, this property would lead to significant 
variations in the lineshape with if it survived in the 
real electronic structure. Our results suggest that this is 
possible only at small Zt . We discuss the band structure 
effects in the following section. 

The panels (B) and (D) show the same as (A) and 
(C) but for smaller wavevcctor k^^ = 0.6 A~^. Compar- 
ison between the right and left thus demonstrates the 
strong dependence of the substrate Green's function on 
the wavevector itself, not just the product kuR\\- The 
two most significant features are that with decreasing 
fco,.: (1) the frequency and damping of the oscillations 
with decrease and (2) the dependence on Zt weak- 
ens. In comparing the two different energies, we assumed 
that the damping constant A (i.e. the tunneling barrier) 
is identical in the two cases. This would be the case 
in metals with identical work functions for states at the 
Fermi energy, in one of which the bottom of the band 
were closer to the Fermi level (smaller k^^). We note that 
kui = 1.2 A^^ corresponds to energy cj = 5 eV in the 
middle of the parabolic band with our parameterization. 
Therefore, the value of A{Rt = 0) = at this energy but 
is negative for smaller energies, e.g. for fc^^ = 0.6 A~-^, 
since in this case there are more high energy continuum 
states repelling the discrete state downward than low en- 
ergy states pushing it up. We also see that the value of 
A at = can change sign with Zt depending on the 
energy u). We note that, since A enters the expression 
for q, the Fano parameter could also be negative and the 
asymmetry of the resonance lineshape could be reversed. 

G. Electronic structure effects and the surface states 
on (111) noble metals 

In the previous section we introduced a simple model 
of G+ based on the unperturbed jellium surface. In gen- 
eral, more realistic behavior of G+ can be obtained from 
electronic structure calculations. Here we discuss qual- 
itatively the electronic structure effects with special at- 
tention to the (111) surfaces of noble metals frequently 



used in STM studies. 

It is well known that (111) surfaces of noble metals 
contain Shockley surface states inside the projected two 
dimensional band gap that forms on these surfaces. 
Both the surface state and bulk wavefunction are given by 
the same general expression (2.11) outside of the metal 
surface. However their overall degree of localization at 
the surface is determined by the position of Cgs, the sur- 
face state eigenvalue, with respect to the band gap edges. 
All other things being equal, the most localized surface 
state occurs when ess is at midgap. As Css moves to- 
wards either band edge, the extension of the evanescent 
oscillatory tail of the surface state wave function into 
the bulk increases, ultimately becoming identical to a 
periodic Bloch function when Css hits the band edge. 
From elementary normalization considerations, surface 
state extension into the bulk and amplitude at the sur- 
face, as reflected in the scale factor (or normalization con- 
stant) for the surface state tails (3.29) extending into vac- 
uum, are intimately related; greater population within 
the bulk means lesser in the surface region. ^^'^^ This sur- 
face state delocalization into the bulk allows for the local 
density of bulk states at the surface to greatly exceed 
that of the surface states, in which case the relative im- 
portance of the surface state in the tunneling current will 
be small near the surface. However, its importance in- 
creases with increasing distance from the surface because 
the bulk states with shorter wavefunction tails are elimi- 
nated from the tunneling. The surface state accounts for 
about 50% of the total signal in typical STM tunnel junc- 
tions in Au(lll)'''^"''^ and is known to be responsible for 
the interference effects observed on these surfaces near 
edges, impurities, and in quantum corrals. It will 
also play a disproportionately important role in the res- 
onance tunneling at large lateral tip-adsorbate distance 
because its contribution to Gq does not decay as quickly 
as that for the bulk states. 

We see the different behavior of the bulk and the sur- 
face states in the STM when we consider the propaga- 
tor Gg for the Shockley state. This is again given by 
Eq. (3.31). However, the k-sum now only extends over 
the 2-D wavevcctor k\^ . Assuming parabolic dispersion 
for the surface state, 7 is given by 

7{Rt,oj) = Jo(fca,i?||)e-^^*'='/' (3.35) 

and 

A{Rt,u,) = ^— depsie)^^^^, (3.36) 

TTps(uj) J„ oj-e 

where, as before, k^^ is the 2-D wavevector of the sub- 
strate state corresponding to energy u. The contribution 
of the Shockley state to the normalized density of con- 
duction states viZt) is given by 

K^t,o;) = e-^^''='. (3.37) 

The propagator Gq for the surface state is essentially 
equal to the Bessel function Jo{ku,R\\) weighted by the 
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exponential e~^^*'^^/^. Therefore the oscihations are not 
damped with increasing Zt and only their overall ampli- 
tude is diminished. Since the surface state on the noble 
metal surfaces (111) have a short kp ~ 0.15 — 0.2 A~^, its 
propagator will have a much longer spatial extent than 
that of the bulk states. The corresponding oscillations 
thus have a spatial period of about 10 times that of the 
Bessel function Jq in Fig. 3(A) in agreement with the ex- 
perimental observation of Friedel oscillations. The contri- 
bution of the surface state to the total current can carry 
information about an impurity on the surface over a long 
distance. The perturbation of the surface state by the 
impurity should persist over several tens of angstroms. 

It is also known^'^ that the spectral weight of the sur- 
face state decreases near surface imperfections. We ex- 
pect the same to be true near the adsorbate. While we 
have explicitly taken into account the interaction of the 
conduction states with the discrete state "a" through the 
adsorbate Green's function G^, all other adsorbate- metal 
interactions, such as potential scattering of the conduc- 
tion electrons from the adsorbate and hybridization of 
the outer shell adsorbate electronic states with the con- 
duction electrons, are neglected in our model. In princi- 
ple, these "residual" adsorbate-metal interactions can be 
included by modifying G"*" and tkp- Although a realistic 
calculation of the system electronic structure is necessary 
to see the effect of the adsorbate on the behavior of G+ 
around the adsorbate, we believe that it will not produce 
oscillatory behavior in G+. In a typical metal, several 
bands with anisotropic dispersion relations eg contribute 
to G+ giving rise to more complicated behavior with no 
single frequency. This will further reduce any oscillatory 
behavior seen in Fig. 3. At larger distance from the ad- 
sorbate, the band structure of the clean surface will be 
reestablished and, as a result, tunneling into the surface 
state. 

Since the importance of the direct tunneling into the 
tightly bound impurity orbital "a" relative to the tun- 
neling into the metal should be weak and decreases with 
increasing iJy, it is useful to study the asymptotic be- 
havior of the conductance Q in the limit to = 0. This is 
equivalent to replacing the Fano parameter q{Rt, Ro;u)) 
by A{Rt; iv) inside Y{Rt, Ro; to), (3.25), in the expression 
for conductance. It then follows that, if the oscillations 
in G+ persist, the lineshape should change with i?|| and 
antiresonances should form at positions where A^ > 7^- 
Using G+ = 7rpse"^/^(A - ij) and 

Im|G+(i?t,i?o)GfG+(i?o,i?t)} = (3.38) 
-7r2p2e-22*/^ { (^2 - A2) Im G^ + 2A7 Re Gf } 

we can write AQeq = Geq — Go at zero temperature by 
replacing {-df{Lo)/du}) by S{lo — V) and using tc{Z) = 
toe~^/^ as 

AGeqiV) = -'^-^tlpMV^ ^ (3.39) 



X lm[G+{Rt,Ro;V) G^{V) G+{Ro,Rt;V)] . 

We see that the resonance in the conductance is a re- 
sult of an interference between different conduction states 
scattering resonantly from the impurity. Its long range 
behavior on the (111) noble metal surfaces is controlled 
by the surface states. Whether the resonance can be ob- 
served at the large distances (> 20 A) depends on the 
spectral weight of the surface state and on its hybridiza- 
tion (~ Vjj2) with the impurity orbital "a" (usually d or 
/). Interesting spatial effects may be realized in system 
with suitable boundary conditions. We believe that Ei- 
gler's quantum mirage of the Kondo resonance inside the 
elliptical corral falls into this category.^** Based on the re- 
sults of the previous section, we do not expect "Friedel" 
oscillations at smaller distances and with period of a few 
angstroms characteristic of the bulk k-vector, although 
we cannot completely rule these out for small Zf. 

IV. DISCUSSION AND EXAMPLES 

The equations (3.24) and (3.27) were derived under 
rather general assumptions. They are suitable as a start- 
ing point for numerical investigations given the necessary 
input from electronic structure calculations. In the rest 
of the paper, we discuss the implications of our theory 
for several specific cases of interest. In all of these cases 
we use our simple model for G"*" based on the jellium sur- 
face and the DOS given by /Os(w) = /Ot(w) = 1 — ui^/D^ 
with D = 5 eV the band half width. In order to elimi- 
nate the exponential fall-off in the tunneling conductance 
with the tip-surface separation and the background dis- 
tortions, we plot the normalized change in conductance 
due to the additional impurity defined as 

AGegiV) = hiGegiV) - go{V))/(27re^pt{Q)tl), (4.1) 

where Qeq is given by (3.27) and Go by (3.28). This is 
equivalent to replacing Y by AY = Y — u in the expres- 
sion for Q. 

Although we were motivated by the experimental ob- 
servation of the Kondo resonance^^'^^ and this work is 
mostly applied to the tunneling through the Kondo res- 
onance, we discuss many of the tunneling properties on 
the simple noninteracting model. We do this primar- 
ily because most of the STM observable characteristics 
of the Fano resonance are common to the single particle 
and Kondo resonances, despite the difference in processes 
that give rise to the two resonances. We wish to point 
out these general features on a model that is conceptu- 
ally far simpler and more familiar to the surface science 
community than the Kondo model, and emphasize that 
the resonances can also be observed in systems with non- 
magnetic impurities with a tightly bound orbital near the 
Fermi level. Finally, the connection with Fano result and 
the consequences of the spatial resolution of the STM 
become more transparent when the same noninteracting 
Anderson Hamiltonian is used. 
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A. Noninteracting adsorbate 

We begin our discussion with an adsorbate-metal sys- 
tem described by the non-interacting Anderson model 
(U = 0). The impurity resonance is characterized by 
its energy eo and the width Tas- The retarded Green's 
function G^{u) for the adsorbate state is 



(a;-eo-ReE„ + ir„^/2)- 



(4.2) 



where Tas{Zo,u}) is 



(3.16) and 

^ is the real 



defined in Eq 

Re E„(Zo,c^) = PEfel^afe(^o)P(^ - 
part of the self-energy for the noninteracting Anderson 
model (not to be confused with the real part of the sub- 
strate Green's function A given in Eq. (3.18)). Following 
Fano we now define the dimensionless energy parameter 
e{Ro,u;) by 



e = 



2{tj - eo - Re E„) 



(4.3) 



We neglect all nonequilibrium effects since they are 
likely to be insignificant for the noninteracting system 
under most experimentally realizable conditions. The 
differential conductance, in lowest order in tc and ta, is 
given by Eq. (3.27) where Y for the noninteracting sys- 
tem takes the form 



Y = Yn = u 



TT^2 • 



(4.4) 



All terms are evaluated at energy w and at the appro- 
priate tip position. We note that lo(0,0,a;) = loo char- 
acterizing the unphysical case of the STM tip in contact 
with surface at the position of the adsorbate (embedded 
in the surface) has the analytic form obtained by Fano 



^00 = 



(g + e)^ 



c2 ' 



(4.5) 



although the inherent energy dependence of q (through 
A(a;)) could distort the pure Fano character of the line- 
shape, even for this "almost atomic physics" STM exam- 
ple. At all other tip positions, the shape of the resonance 
will be described by the more general expression (4.4), 
which is equivalent to a sum of a "Lorentzian" plus a 
"Fano" profile, 



Ycs = 



1 + e^ 



(g^ + e)^ 
' 1-Ke2 



(4.6) 



with v' = v + q^ —^"^ — (^J^ jv and g' = q^jv. Decompo- 
sition of Yq in this manner may be useful when analyzing 
experimental lineshapes. Note that the additional intrin- 
sic energy dependences of v and 7 could further distort 
the standard lineshape. 
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FIG. 4. Differential conductance AQeq as a function of the 
strength of the direct tunneUng matrix ta (in units of tc). 
The model parameterization is given in text. The two curves 
correspond to two different impurity level energies eo = 5 eV 
(solid) and 2 eV (dashed) from the bottom of a symmetric 
band 10 eV wide. The two insets show the model density 
of states ps for the conduction electrons (upper panel) and 
the imaginary (soUd) and real (dashed) part of Ga(tt') (lower 
panel). 



1. Lineshape dependence on electronic structure and on the 
relative strength of ta and tc 

We first show (Fig. 4) the dependence of AQ^q on the 
ratio {ta/tc) for .R* = and a resonance at the Fermi 
level. The solid line corresponds to a resonance at the 
center of a parabolic band (symmetric around its center), 
i.e. eo = 5 eV from the bottom of the band and the 
dashed line corresponds to a resonance at eo = 2 eV from 
the bottom of the band. The two energies correspond 
to the band energy Cfc with wavevector k = 1.2 and 
A; = 0.6 A-i in our jellium model (Fig. 3). The resonance 
width is r^jg — 0.2 eV in both cases. At zero temperature, 
from Eq. (3.27) and (4.4) we write 



Qeq{V) 



2e^ 
~h 



Ps{V)T,pt{Q)tlYo{V) 



(4.7) 



In order to make connection with the Fano result, we 
plot the conductance for small tip-metal separation with 
the tip above the adsorbate {Rt = 0). The lineshape 
Yq is then given by the Fano formula (4.5). The plotted 
quantity A.Qf,q{V) in Fig. 4 is then given by 



AGeqiV) = Ps{V) 



1 



(4.8) 



The Fano parameter q depends not only on the ratio 
(ta/tc) but also on energy and electronic structure. We 
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see this most clearly in the first panel where ta = 0. The 
resonance placed at the center of the band produces a 
symmetric dip in AQ characteristic of g = 0, whereas the 
resonance at eo = 2 eV has an asymmetric lineshape due 
to the negative contribution from A to g (see Fig. 3). Its 
lineshape actually becomes symmetric at finite value of 
ta- The value of ta inside each panel is given in units of tc- 
The inset in the upper panel shows the model density of 
conduction states ps and the lower panel inset shows the 
spectral function pa = — ^ImG^ (solid) and ReG^/w 
(dotted) for the level at the center of the band. 

As the strength of the direct tunneling ta increases 
with respect to tc, the resonance develops its characteris- 
tic asymmetric shape and, eventually, at large ta/tc ^ 1 
it acquires the shape nearly indistinguishable from that 
of the impurity spectral function /9o(w). With increas- 
ing tip-adsorbate separation, the signal from the reso- 
nance must disappear as both the tunneling element ta 
and G~^{Rt,i^) tend to zero. The differential conductance 
is then determined by the density of states of the clean 
surface. This property is not present in the Fano expres- 
sion. We now discuss this behavior. 



2. Lineshape dependence on the tip-surface separation 

Using the same model system as in the previous section 
with the resonance at the center of the band (eo = 5 eV 
from the bottom of the band), we demonstrate the de- 
pendence on Zt (with = 0) in Fig. 5. We make the fol- 
lowing model for the tunneling matrix element ta{Rt,Ro) 
and tc{Rt)- The exponential fall-off of the metal and 
adsorbate wavefunctions is controlled by different decay 
constants. The adsorbate state Va is tightly bound es- 
pecially for narrow resonances of interest here. The con- 
duction electron wavefunctions, on the other hand, typi- 
cally belong to the outer s or p orbitals and have longer 
tails into the vacimm. As a consequence, the ratio ta/t^ 
and thus also the Fano parameter q, changes with Zf. 
In order to incorporate this property, we use the matrix 
elements (2.17) and tc = t^e'^^l^, where a = 0.75 A, 
A « 0.9 A, io = 25 meV, and ta = O.ltc at = 2 A. 
Under these conditions, the q parameter tends to zero 
with increasing Zf. 

The panel (a) shows the normalized AQcq for this 
model. The lineshape undergoes only moderate changes 
with Zt within the experimentally relevant range. We ex- 
pect this to be a general property. In order to understand 
the behavior, we discuss the lineshape dependence on tip- 
surface separation conceptually in terms of two contribu- 
tions: (1) different decay constants for the discrete ipa 
and metal ipk states at the Fermi level, and (2) different 
decay constant for metal states at eps with different fcy. 
We separate the observable consequences of these two 
effects in panels (b) and (c). The first contribution pro- 
duces changes in q due to the changing relative strength 
between ta and tc- We demonstrate this in Fig. 5(b) 



where only this contribution is taken into account by set- 
ting Zt = inside the substrate Green's function , or 
equivalently by setting 7 = z/ = 1 and A = in Eq. (4.4). 
This limit does not correspond to a real situation and 
does not lead to the correct Zt 00 limit. It is shown 
here merely as an example of the contribution (1) to the 
Zt dependence of the tunneling conductance. With our 
parameterization, this case is identical with q changing 
from q « 0.8 at = 3 A to g « 0.2 at Zf = 9 A. As 
Zt increases further, g — > and the resonance becomes 
symmetric. However, we see that the normalized conduc- 
tance AQeq does not vanish in the limit Z^ — > 00. 
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FIG. 5. Differential conductance AQeq vs. Zt for a model 
described in the text, (a) the complete dependence on the tip 

distance Zt\ (b) assumes Zt = Q inside and thus neglects 
the wavevector dependent effects; (c) is same as (a) but with 
ta = 0, i.e. it only includes the effect of decreasing spatial 
sensitivity incorporated through the substrate Green's func- 
tion . The vertical scale is arbitrary but identical in all 
panels. 

Fig 5(c) takes the second contribution (2) into account 
while leaving out the first one. We chose ta/tc = = 
constant, which would be the case if there were no di- 
rect tunneling into the discrete state. In this extreme 
limit, any changes in lineshape arc a consequence of the 
varying weight that different fcy metal states play in the 
tunneling at different Zt- This occurs because the k\\ = 
metal wavefunctions given by Eq. (3.30) have the great- 
est extension into the vacuum and as a result the spatial 
resolution of the tip decreases. Therefore the signature 
of the resonance in Q decreases even after normalization 
of the current for different Zt as the ratio AY/v — > 
with Zt ^ 00. Fig 5(a) shows the combined effect of the 
contribution and represents realistic conditions. It ac- 
counts correctly for the changing lineshape, as well as its 
disappearance. We again emphasize that realistic band 
structure is desirable for making quantitative statements. 
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Finally, wc comment on the experimental issues. It is 
clear that the lincshapc dependence on Zt will be observ- 
able only if it can be studied over a reasonably large range 
of Zf, this being limited by the experimental resolution 
and detection capabilities. The most favorable case is 
one in which the direct tunneling ta into the resonance 
is strong at small Zt, i.e. q is large, and the difference 
in relevant decay constants for the adsorbate and metal 
electrons, a and A respectively, is large. This is not the 
case in the experiments^^'^^ where q is small. Therefore, 
we do not expect significant changes in the lineshape with 
Zt in these systems. Since our model is based on realistic 
parameterization, we expect the behavior shown in Fig. 5 
to serve as a guide for order of magnitude estimates for 
the spectral dependence on Zt. 

The direct effect of the STM tip on the system and 
thus also on the lincshapcs is not taken into account here. 
This issue is discussed in section IV B 3. 



3. Lineshape dependence on the lateral tip position 

As wc already discussed in section IIIF, the resonance 
lineshape depends on the relation between the spatial 
dependence of the direct tunneling and the propagation 
of the adsorbate- induced perturbation through the metal. 
Fig. 3 shows that the direct tunneling into the resonance 
(ta) is expected to fall off faster than the perturbation. 
Therefore, at large , the lineshape will be given by the 
ta = limit conductance. We show the dependence of 
AQeq on in Fig. 6. We do this again for the model 
described in the previous section (Fig. 5) with ta = 
and a resonance at the center of the band, i.e. eg = 5 eV 
from the bottom of the band which corresponds to k^^ = 
1.2 in Fig. 3. The solid line corresponds to = 5 A 
and the dashed line to Zt = A. 

The unphysical case oi Zi = Q (dashed) is shown to 
emphasize the possible consequences of the oscillations 
in Gq displayed in Fig. 3. We chose the lateral tip posi- 
tions in the figure to coincide with thc^ nodes and zeros of 
A and 7 to show the dramatic changes in the lineshape 
with due to the oscillations in A and 7. Since the 
spatial decay of the oscillations is small at Zt = 0, the 
sequence of resonances and antiresonances appear in the 
range i?|| S (0, 10) A. The possibility for such antireso- 
nances is discussed implicitly in the work of Kawasaka^^ 
and explicitly by Schiller and Hershfield^^ . However, this 
behavior is not observed in the experiments by Madhavan 
et al.'^^ and Li et al.'^'^ due to the smoothing of the elec- 
tronic structure with increasing distance from the surface 
that we discussed in IIIF. 

In fact wc would not expect the dramatic variations in 
line shape with reported by Schiller and Hershfield'^^ 
to be observed. The reason is apparent from the behavior 
of G+ as a function of Zt (Fig. 3). As the tip distance 
from the surface increases the oscillations axe destroyed 
by the increasing weight of the lower frequency (small 



) components at larger Zt interfering destructively with 
those given by fc^. For this value of k^, the oscillations 
are effectively damped when > 5 A and the shape of 
the resonance does not change significantly as shown by 
the bold line in Fig. 6. We expect that band structure 
effects will suppress the oscillations even further. 
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FIG. 6. The normalized differential conductance AQeq as 
a function of the lateral tip position for the same model 
as in Fig. 5: Zt = 5 A(solid) and Zt = A(dotted). The 
contribution of the surface states to is not included here. 

We also find that the spatial extent of the resonance in 
the spectrum should decrease as the STM is retracted, as 
long as the signal is due to the bulk states. At Zj = A, 
the resonance is still visible at R\\ ~ 10 A but only to 
about 4 A at Zf = 5 A. This is a somewhat shorter 
distance than that found experimentally for Co/^u(lll) 
and Ce/ylg(lll)^^'^^. Although the Fermi wavevector 
kp ~ 1-2 A~^ used in Fig. 6 is close to the free electron 
value of kp for the noble metals, the disagreement is not 
surprising since we made no real attempt at realistic elec- 
tronic structure description. Smaller values of kp would 
increase the spatial extent as would smaller values of Zt 
and A. 

Interestingly, the Zt = [dotted] lineshape progres- 
sion shown in Fig. 6 is qualitatively similar to the family 
of lineshapes that would be expected from surface state 
propagation, but with i?,|| , the lateral tip-adatom sepa- 
ration rcscalcd upward by nearly an order of magnitude. 
This claim is based on the qualitative similarity between 
the bulk Gq at Zt = and the surface Gq . The bulk 
7(i?l|,eF) = jo{kpR\\) [the dashed curve in Fig. 3] at 
Zt = the analogous surface state 7(i?||, e^) = Jo{kpR\\) 
[dotted curve in Fig. 3] both exhibit long range oscilla- 
tions unlike the bulk state at Zt > SA. However since 
kp ~ (0.1 — 0.2) A^-'^ for the surface state band, Jo{kpR\\) 
shown in Fig. 3(A) for kp = 1.2A~^ should be plot- 
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ted with this sniaUer kp when referring to actual noble 
metal surface state bands, in which case the observable 
i?ll -dependent lineshape evolution in Fig. 6 would still 
be representative, but with rescaled by the factor 
1.2/0.15 = 8. From this it is easy to appreciate that 
the dramatic lineshape variations will occur mainly at 
very large lateral separations. Clearly, realistic electronic 
structure calculations are necessary to answer the more 
quantitative questions. 

At large values of and for broader resonances, an 
additional mechanism for distortion of the lineshape is 
possible if the relative change in the length of in the 
energy range given by the resonance width (F) near eo 
is large. In principle, this gives rise to the possibility of 
the oscillatory behavior of Jo with respect to energy ex- 
hibiting itself in the spectrum. However, the conditions 
for this effect in G(lu) would require the relative change 
Ak,,/ke « + F/2eo - a/1 - r/2eo to be > 27r/fc,i?||. 
In this case, the argument of the Bessel function in (3.32) 
will vary over several periods starting at a small value at 
the bottom of the resonance. It is however clear that the 
width of the resonance and the distance from the impu- 
rity would have to be much larger than in the recent ex- 
periments with Kondo impurities. Moreover, these oscil- 
lations could only be observed through the surface state 
due to the damping of oscillations in the bulk. 
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FIG. 7. A^eq as a function of temperature for Zt = 5 A 
and the same system as in Fig 5. Temperature and energy 
are given in units of the resonance width, V. 

4- Temperature dependence of differential conductance 

The temperature dependence of the tunneling conduc- 
tance through a narrow resonance is an important is- 
sue especially in considerations of tunneling through the 
Kondo resonance which itself is temperature dependent. 
We revisit this issue in the section IV B 2 on Kondo ef- 



fect. Here, we demonstrate the effect of the broadening 
in the Fermi function on the spectrum when the impurity 
spectral function is independent of temperature. 

We assume the density of tip states and the substrate 
conduction electrons to be constant on the scale of the 
width F near the tip Fermi level and around the reso- 
nance. The temperature dependence of the conductance 
then follows from equation (3.27) with (4.4). We consider 
only the case of the tip directly above the adsorbate for 
simplicity and write 

g{v) = '^-^^p,p,p^d^{-^-Mp.^ t%{R„uj). 

(4.9) 

It follows from this expression that the differential con- 
ductance only depends on the temperature of the STM 
tip. This is a consequence of the assumption that the 
density of tip states is constant in the relevant energy 
range and the tunneling barrier is not modified by the 
bias. The differential conductance is then a function of 
the spectral density of the substrate states and is inde- 
pendent of their occupation. Generally, if the DOS in 
the STM tip varies significantly on energies ~ Fq^, the 
substrate temperature would also enter. If the density of 
states in both the substrate and the tip were constant, 
no temperature dependence would be observed. It follows 
from Eq. (4.9) that the temperature T must be ~ F, in 
order to have a significant effect on the conductance. 

We show the temperature dependent Qgq in Fig. 7. The 
spectral density ^0(0;) with the narrow resonance at the 
Fermi level (the same system parameters as in Fig. 4) 
is independent of T. The STM tip is directly above the 
resonance with Zf = 5 A. In the range T < 0.1 F, no 
temperature dependence is noticeable. However, when 
T ~ (0.2 — 0.5) Fas the differential conductance shows a 
rather strong dependence on temperature, and at T F 
the sensitivity of the STM to the resonance disappears. 
The temperature dependence in Fig. 7 comes entirely 
from the broadening of the Fermi function of the tip. 

B. Tunneling into Kondo resonance 

In the previous section IV A, we discussed the STM 
conductance in tunneling through a nonintcracting im- 
purity [[/ = in Hamiltonian (2.1)], frequently referred 
to as the resonant level model (RLM). We now turn to 
the case of magnetic impurities and tunneling through a 
Kondo resonance. We begin with the case of a weak 
tip)-metal coupling. However, for the; Kondo systems 
this assumption is more restrictive than for the RLM 
model, and for this reason, we later take advantage of 
our noncquilibrium approach to account for the direct 
effect of the tip on the impurity spectral density, while 
still neglecting the tip's effect on the metal states. 
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1. Conceptual and theoretical approach 

Since our earlier derivation of the current and con- 
ductance is valid for arbitrary interaction [?7 in 
Eq. (2.1)], the final results (3.24) and (3.27) also hold 
in the Kondo and mixed- valent regimes of the Anderson 
model (i.e. U ^ A). The properties of the adsorbate 
enter through the Green's function Ga- The problem is 
thus reduced to finding the one electron Green's function 

Ga. 

However, we first consider the tunneling for a spin 1/2 
(a = cr) impurity in the Kondo limit, {eps — eo) T and 
(eo — €fs + U) » r. The Kondo resonance has a very 
small weight and is due to spin fluctuations. The possible 
tunneling channels in this case are shown in Fig. 8 as 
processes (1) and (3). The system can be described by 
the Kondo Hamiltonian in this limit 

Hs{Zo) = J2 efecLcfe, + ^P^^cpa + (4.10) 

kcr V<r 
+ I^{*fep(^t)cLcpa+H.C.} + 

kk' uu' 

)-S + H.c] 

kpaa' 

where the first three terms were also present in the to- 
tal Hamiltonian introduced in section H and describe the 
unperturbed metal and tip states and the coupling be- 
tween the two. The remaining terms give rise to spin 
fluctuations in the presence of the magnetic impurity. 
The terms with couplings Js and Jt correspond to the 
exchange interaction of the local spin with the substrate 
and tip electrons, respectively. The last term (Jst) corre- 
sponds to the effective tip-substrate exchange interaction 
in which charge is transported between the tip and the 
surface. This Hamiltonian can be obtained from Htot of 
section H using the Schrieffer-Wolf transformation which 
relates Jg, Jt and Jst to Va and ta- For the symmet- 
ric Anderson model, Jg = 4:Vq /U, Jgt = ^taVa/U, and 
Jt = ^t\/U. Using the continuity equation (3.1), the 
current is 

/ = y MjI^'^P^'^Lcpa) + Jst Y ^Aa^aa'Cp^' ' S) + 
kpcr kpaa' 

+ Jt J2 (4^WCpV' • S)}. (4.11) 

pp' aa' 

The first term is identical with the first term in 
Eq. (3.3). In the lowest order of the tip-system couplings 
{tkpT JstT Jt)i the third term does not contribute. The 
first term corresponds to the direct tip-substrate tunnel- 
ing channel - process (1) in Fig. 8 - which includes the 



scattering of conduction electrons from the local moment. 
The second term corresponds to the direct tunneling into 
the magnetic impurity - process (3) in Fig. 8. We note 
that the spin flip scattering that gives rise to the Kondo 
effect is a higher order process. In the lowest order, the 
channel (1) and the spin-flip component of (3) do not give 
rise to interference because the final states have different 
spin states. The lowest spin-flip process that does inter- 
fere with (1) is of second order in J and proportional to 

Js^st' 

In the limit of large tip-metal separation, equivalent 

to the condition (Js ^ Jst ^ Jt), the third term in 
Eq. (4.11) as well as higher order contributions from Jgt 
are neglected and all other exchange processes are in- 
cluded in principle. This is equivalent to assuming that 
the state of the metal-adsorbate system is determined 
only by Js and is unaffected by the presence of the tip. 
Theoretically, the problem then reduces to flnding the 
spectral properties of the system without the tip and us- 
ing them in the expansion for tunneling via the two terms 
in Eq. (4.11). 

As the system parameters move away from the Kondo 
limit - that is either eo shifts towards eps or U becomes 
smaller - valence fluctuations appear. The Kondo reso- 
nance is then due to both the spin and charge fluctua- 
tions. The separate energy scale due to the spin fluctu- 
ations eventually disappears in the mixed-valent regime 
and the Kondo peak merges with the broad resonance 
centered at cq. In the intermediate regime, where both 
charge and spin fluctuations coexist on the impurity, an- 
other tunneling channel exists. This channel is denoted 
by (2) in Fig. 8. It also includes the contribution from 
higher order non-flip processes similar to (3). We study 
the system in this regime with the Hamiltonian defined 
in section II. 




FIG. 8. Possible scattering channels for an electron tunnel- 
ing from tip to metal through a magnetic impurity adsorbed 
on the surface. 

We adopt the slave-boson technique of Goleman^^ 
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and find the adsorbate Green's function using the non- 
crossing approximation (NCA)^^. Following the theory of 
section III, the final expression for current in Eq. (3.5) is 
valid, as well as all the consequent steps and approxima- 
tion in III. We insert the solution for the Green's function 
Go of the ([/ = oo) interacting system in Eq. (3.27). This 
is equivalent to including the three tunneling channels in 
Fig. 8 to lowest order in the tip-system couplings. We 
now turn to the discussion of the results based on this 
approach. 



2. Results for large tip-substrate separation 

In order to model Co/^u(lll) studied both experi- 
mentally^^ and theoretically^^, we choose a parameter- 
ization that gives the Kondo temperature Tk « 70 K 
appropriate for the system. Our simplified model has de- 
generacy N — 2 with no orbital degeneracy, band width 
2D = 10 eV, and the adsorbate level at eo = 0.75 eV 
with the width 2r = 1 eV (the width of a multiplet with 
an occupied level is AT rather than F!). We show the 
corresponding spectral function and the real part of Ga 
in the inset of Fig. 9. 
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FIG. 9. The spectral lineshape for a model Kondo system 
described in the text over large bias range that includes tun- 
neling into the broad resonance at 0.75 eV below the Fermi 
level. Each panel corresponds to different q{eFs) at the Fermi 
level. 

Fig. 9 shows the spectral properties of the system using 

AQ^g over the whole energy range of the conduction band 
and for diflFerent values of q at the Fermi level. The spec- 
trum contains information about the broad resonance at 
eo = 0.75 eV below the Fermi level, as well as the promi- 
nent feature due to the Kondo resonance at zero bias. We 
show the large bias voltage results only for completeness 
since we do not expect the STM experiments to be able 



to provide spectroscopic information about the system 
over the whole energy range shown. 

The resonance lineshapes both in Co/Au{lll) and 
Ce/Ag {111) correspond to small values of q. Madhavan 
et al?^ fitted the observed resonances to Fano lineshapes 
with q ~ 0.7. Our best fit would give approximately 
the same value of q. In the case of Ce/Ag{lll), the ob- 
served feature is an almost symmetric antiresonance cor- 
responding to q ^ 0. Due to the contribution from the 
substrate electronic structure to g, its value cannot be 
directly used to make quantitative statements about the 
relative strength of the tunneling into the discrete state 
d (/) with respect to that into the continuum. However, 
in agreement with Li et al?^ and Lang^^ we conclude 
that the STM probes mostly the sp wave functions and 
the tunneling into the /-orbital is rather weak at the tip- 
adsorbate distances used in the Ce/Ag{lll) experiment. 
The resonance is mostly the result of interference between 
conduction electrons scattering from the impurity. The 
larger value of q in Co/Au{lll) indicates stronger con- 
tribution from the coupling of the STM to the d orbital. 
This is expected because the 3d orbital is not as tightly 
bound. 

The recent work of Kawasaka et al.^^^^^ deals with the 
spatial and spectroscopic profiles of the Kondo resonance. 
They begin with the Tersoff-Hamann^^ expression for 
current (3.12) and use the local density of states given 
by (3.14). They insert the self energy correction in the 
Green's function Ga due to the intra-adsorbate Coulomb 
correlations using perturbation theory (T > Tk) and Ya- 
mada's expansion in [/ (T < Tk) to study the tem- 
perature dependence in the whole temperature range. 
They neglect the additional temperature effects due to 
the Fermi surface broadening, replace {—df/duj) by the 
delta function, and evaluate the conductance at the tip 
bias. 

One of the main conclusions of their work^^ is that the 
calculated temperature dependence of the resonance in 
the differential conductance is indicative of the temper- 
ature dependence of the Kondo resonance itself. They 
show results at the experimentally relevant low temper- 
atures for Co/Au{lll) and Ce/ Ag{lll) in the range of 
temperatures (T < O.ITk)- They find a rather weak tem- 
perature dependence, due entirely to the temperature de- 
pendence of the spectral function pa- It is qualitatively 
the same and comparable in magnitude with that found 
in Fig. 7 for a temperature independent resonance of the 
noninteracting system for temperatures T < F. In our 
case, the temperature dependence in /S.Qeq is the conse- 
quence of the Fermi surface broadening in the STM tip. 
Therefore a careful deconvolution is necessary even at 
these low temperatures to extract information about the 
temperature dependence of the Kondo resonance. The 
other possibility is to eliminate variations in the Fermi 
surface broadening of the tip. 

We show the temperature dependence for a Kondo sys- 
tem in Fig. 10. Since the validity of our approximation 
is limited to temperatures of order Tk and higher, we 
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show our results only in this temperature range. Panel 
(c) shows the temperature dependence one would observe 
with the tip at T = K and with varying substrate tem- 
perature, i.e. when only the temperature dependence 
of the spectral function is taken into account. Panel 
(b) assumes the substrate is at a constant temperature 
T = Tk, which determines the shape of the Kondo reso- 
nance, while the tip temperature is varied. We see that 
the two contributions produce a very similar broadening 
of the Fano resonance. Only a close look can uncover 
the difference. Panel (a) shows the combined effect when 
the tip and substrate are kept at a common temperature. 
Obviously, it would be difficult to determine the contri- 
bution from the broadening of the Kondo resonance. 
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FIG. 10. Temperature dependence of the tunneling con- 
ductance through the Kondo resonance, (a) the total de- 
pendence includes the temperature dependence in the spec- 
tral function and the Fermi surface broadening; (b) tempera- 
ture dependence duo to the Fermi level broadening and with 
spectral function given by its T « Tk value for all tempera- 
tures; (c) temperature dependence of only the spectral func- 
tion without the Fermi surface broadening (replaced by delta 
function) . 

In addition to the temperature effects just discussed, 
Kawasaka et al.^^'^^ also predicted the existence of weak, 
long range oscillations in the current as a function of 
the lateral tip position. The particular long wavelength, 
long range character of these predicted oscillations are 
a consequence of their assumption that the tip-to-mctal 
tunneling takes place into the surface states of the (111) 
noble metal surfaces. The observed resonances^^'^^ are at 
variance with these expectations. On the other hand, the 
limited spatial extent of the resonance observed at lat- 
eral tip positions up to 10 A is consistent with the rapid 
spatial decay determined by the bulk Gq. No significant 
changes in the resonance lineshape are expected on this 
length scale (sec discussion in section IIIG and Fig. 6). 
The surface state would, however, be responsible for line- 



shapes variations on larger length scales of order 20 A. 
The fact that no resonance is observed at such a distance 
from the impurity indicates that the surface state con- 
tribution is indeed weak. Since Friedel oscillations have 
been observed over long tip-impurity separations, we be- 
lieve that a weak tunneling resonance most likely persists 
in the conductance over comparable distances, but more 
sensitive experiments are necessary. In this case, changes 
in the lineshape with are expected. However, unlike 
Schiller and Hershfield^^, we do not expect variations in 
the lineshape due to the dominant contribution from the 
bulk states on the length scale of < 5 A as discussed in 
section IV A 3. 



3. Nonequilibrium and hybridization effects at small 
tip-substrate separation 

In typical STM experiments, the tip-substrate separa- 
tion can be varied from the point of contact where the 
tunneling resistance i? is a few lOOkQ to distances where 
R ~ IGfl. Experimental constraints limit the STM use- 
fulness to the near Fermi level spectroscopy - especially 
at small - because of exponentially increasing tunnel- 
ing currents with bias. However, it is likely to be possible 
to investigate the Kondo resonance which only requires 
biases of the order of ~ 10 meV - with very small tip- 
adsorbate separations. It is therefore useful to analyze 
the physical consequences of the small tip-metal separa- 
tion on the resonance in tunneling conductance. 

In this case, nonequilibrium effects, as well as the tip- 
adsorbate interaction, become important in the spec- 
troscopy of Kondo systems. First of all, as Fat increases 
and becomes a significant fraction of Tas at small dis- 
tances, the tip- adsor bate hybridization will contribute to 
the width F of the resonance and to the renormaliza- 
tion of the level eo- As a result, the Kondo temperature, 
which depends sensitively on F and cq, will change. This 
could be particularly important for systems with very 
low bulk Tk, such as Fe/Au with ~ 1 K. The Kondo 
temperature for an impurity adsorbed on the surface of 
the metal is even lower than its bulk Tk because the 
lower coordination number for the adsorbate makes the 
width F narrower. If Tk <^ T, the Kondo resonance 
will not be observed. In certain systems and in the right 
temperature regime, it may be possible for the Kondo 
resonance to reappear at smaller tip-adsorbate distance 
as a result of the increased hybridization. This could also 
be achieved by incorporating the adsorbate into the top 
surface layer. The recent study of transition- metal im- 
purities at the surface of gold^^ did not find any sign of 
the Kondo effect in V, Cr, Mn, or Fe. We believe that, 
in the case of iron, this is due to the low Tk and may 
be an example of a candidate system for the conditions 
discussed here. On the other hand, the smaller impurity- 
metal hybridization at the surface can lead to magnetic 
behavior for systems which are nonmagnetic in the bulk. 
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such as Ni/Cti. There is a possibiUty for observing the 
transition between magnetic and nonmagnetic behavior 
on a single system induced either by embedding or by the 
proximity of the STM tip. 

We show an example of the changing Tk with hy- 
bridization in Fig. 11, where the spectral function Pa{'^) 
is plotted at zero bias as a function of the partial width 
Tat, i-e. tip- metal separation for a model system. We 
choose D — b eV, Vas = 0.25 eV, = — 1 eV, and 
T = 30 K. The Kondo temperature for this model in the 
limit ta = 0, is ~ 30 mK, much smaller than the 
temperature T. Therefore the Kondo resonance in the 
spectral function is very weak. When the tip is brought 
closer to the adsorbate, the Kondo resonance acquires 
more spectral weight as the Kondo temperature increases 
to Tk ~ 100 mK at Tat = O.OlFa,, Tk 0.2 K at 
Tat = OMTas, and Tk ~ 1.5 K at Tat = 0.25Tas. Based 
on the justifications in Appendix A, we neglected the ef- 
fect of the direct metal-tip interaction on the spectral 
function and treat the effect of the tip as another hy- 
bridization channel for the impurity state. 
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FIG. 11. The spectral function paii^) at zero bias as a func- 
tion of the tip-substrate separation (defined in terms of Fat)- 
Dotted line= occupied density of states, bold line = spectral 
density of states. 

The effect of varying hybridization - due to either the 
presence of the STM tip or due to embedding or chang- 
ing the environment of the adsorbate - on the tunnel- 
ing resonance depends on the relation between Tk and 
T. For instance, when Tk <C T, increased hybridization 
would produce stronger (sharper) tunneling resonance of 
the same width since the spectral weight in the Kondo 
resonance increases while its width remains almost con- 
stant until Tk ~ T. Also the experimental resolution 



is limited by temperature in this regime. When, on the 
other hand Tk > T , additional hybridization would not 
only increase the spectral weight in the Kondo resonance 
but also its width. The two cases should thus be distin- 
guishable experimentally from each other and from the 
possible lineshape variations with Zt as a result of chang- 
ing q. 

The second important effect of the strong tip- adsorbate 
interaction is the breakdown of equilibrium relations at 
finite bias such as the fluctuation-dissipation theorem 
G^{uj) = f{u))pa{u)) which consequently cannot be used 
in deriving the expression for current (3.24). This is true 
in general because the electron occupation of the tip, 
metal, and adsorbate electrons will no longer be thermal, 
i.e. will not be given by /t(a;) and /s(w) but rather will 
be characterized by a nonequilibrium distribution pro- 
duced by the injected tunnel electrons. The differential 
conductance is no longer proportional to the local den- 
sity of states and cannot be obtained using Eq. (3.27). 
In Kondo systems, the hot electrons not only modify the 
electronic distribution on the impurity, but also modify 
the Kondo resonance itself. 

This is shown Fig. 2 where the spectral function and 
density of occupied states is plotted at selected bias volt- 
ages in the limit of \tpk\ <C \tpa\- The impurity has a 
resonance at eo = — 1 eV below the Fermi level and total 
width r = 0.5 eV produced by the hybridization with 
both the tip and the substrate with the partial widths 
Tat = O.lFos- Temperature is of the order of Tk in this 
example. We see that the Kondo resonance broadens 
even more with increasing bias. This is due to the in- 
crease in the rate of incoherent scattering by ~ eVa/T 
- an effect similar to temperature. At the same time, 
the electron occupation develops a non-thermal profile 
due to the large tip-adsorbate current. This is partic- 
ularly visible for negative biases where the density of 
states is larger. Fig. 2(a) shows the equilibrium spec- 
tral function (dotted) and the electron population on the 
resonance (solid bold). The equilibrium spectral density 
is shown (dotted) in all panels. In addition, the spec- 
tral density (solid) and occupation (bold solid) are shown 
for the biases indicated in the figure by the labeled ar- 
row. If the coupling to the tip were comparable with the 
metal-adsorbate hybridization, a double peak structure 
would develop. This has been predicted by Wingrcen and 
]y[gjj-5i,52 ^YiQ context of the nonequilibrium Kondo ef- 
fect in quantum dots, also discussed by Plihal et al.®^. 
We see the onset of the double peak structure in panels 
(b)-(d) where a small cusp develops at the chemical po- 
tential of the tip. In summary, the bias has a significant 
effect on the spectral density even when Tat ~ O.lFag. 

We show the tunneling current (right) and the cor- 
responding differential conductance (left) in Fig. 12 for 
this model of Kondo impurity and for the STM geom- 
etry defined by Tat = O.lFa^ and \tpk\ < \tpa\. The 
panels correspond to q = 0.6, q = 1.2, and q = 2.4, re- 
spectively. The current on the right is calculated using 
Itot = leg + SInon of section III D with 6Inon given by 
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Eq. (3.26). The differential conductance Qtot on the left 
is obtained by differentiating the results displayed on the 
right. It cannot be calculated from the expressions in 
the text since the dependence of /Oa(w) and G^{u)) on 
the bias voltage modifies the contributions to the cur- 
rent in a wide energy range, and Qtot is not related in 
simple terms to the properties at the Fermi level of the 
tip. We compare Itot and Qtot (circles) with leg and Qeq 
(solid). The equilibrium quantities were obtained in the 
lowest order in tap and tkp and with the equilibrium Ga 
of Fig. 2(a). 

We see that the broadening and disappearance of the 
Kondo resonance with increasing bias at strong tip- 
adsorbate coupling is weakened in the nonequilibrium 
calculation of Qtot, because the contribution SInon com- 
pensates partially for the spectral function effect. The 
most consistent effect on the lineshape for various values 
of q is the suppression of the resonance maximum and as 
a consequence a more symmetric appearance. This be- 
havior is qualitatively different from both the hybridiza- 
tion effect and that of the changing Fano parameter q 
- due to different decay constant of the impurity and 
metal states. Although the dependence of the tunneling 
resonance on the tip-substrate separation Zt will contain 
all three contributions, the hybridization and nonequi- 
librium contribution should only be important at ex- 
tremely small tip-adsorbate separations. The variations 
in q should not be important as it depends on the dif- 
ference of the wavefunction tails and the two remaining 
contributions should leave distinguishable signatures in 
the tunneling resonance. It remains to be seen if the 
nonequilibrium condition play an important role in the 
tunneling between STM tip and Kondo impurity. 
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FIG. 12. Differential conductance (left) and current (right) 
for a Kondo system defined in the text and with the spectral 
function displayed in Fig. 2. Solid lino is the calculation in the 
lowest order in tkp and tap and with the equilibrium spectral 
function of Fig. 2. Circles correspond to the nonequilibrium 
Qtot = dItot/dV and Itot for Fat = O.lFas. Each row corre- 
sponds to a given value of q at the Fermi level. 



Finally, we note that the limit \tpk\ ^ \tpa\ discussed 
here in connection with the nonequilibrium effects is not 
appropriate for the recent STM experiments, where tpa is 
likely much weaker than tpk, even though the importance 
of tpa will increase relative to tpk with decreasing Zt. We 
will address the more general case in a future work. 



V. CONCLUSIONS 

We used the Keldysh-Kadanoff method to study the 

spectroscopic features of adsorbate resonances in the 
STM tunneling experiments. The central results of our 
theory are the general expression for current Eq. (3.5) 
to all orders in the tunneling matrix elements and its 
equilibrium limit (3.11). Both are valid for arbitrary 
intra-adsorbate electron correlations and thus apply to 
both noninteracting ([/ = 0), as well as magnetic (large 
U) systems. The discussion of the Fano resonances is 
based on additional approximations for the tunneling 
and hybridization matrix elements that lead to expres- 
sions (3.24) for the tunneling current and Eq. (3.27) for 
differential conductance in the lowest order in the tip-to- 
system tunneling matrix elements ta and tc, i.e. at large 
tip-surface separation, and the nonequilibrium correction 
to the current in Eq. (3.26). 

In the equilibrium limit, our theory of the tunneling 
current and conductance differs from the standard the- 
ories of STM, in that the dependence on LDOS is re- 
placed by a tip specific quantity related to the LDOS 
(Eq. 3.9). The current is expressed entirely in terms of 
the adsorbate Green's function G„ , the tip density of 
states, the tunneling matrix elements and the substrate 
Green's function. We used the formulation to study the 
resonance lineshape as a function of temperature, tip- 
substrate separation, and lateral tip position. We sum- 
marize our findings as follows. 

(1) The role of impurity state resonances in tunneling 
can be discussed in terms of two limiting cases. When 
direct tunneling across the barrier is weak, the resonance 
within the barrier provides an additional tunneling chan- 
nel and can significantly enhance the tunneling current. 
This is the case of quantum dots in Coulomb blockade 
regime. If on the other hand, the tunneling into the con- 
tinuum is strong, the presence of an "impurity" state 
could suppress the tunneling current due to the addi- 
tional scattering of the conduction electrons in the metal 
from the impurity, i.e. increased resistance. The tunnel- 
ing into the Kondo resonance in the recent STM experi- 
ments seems to be closer to the latter limit. 

(2) The information about electron correlations and 
the Kondo resonance enters the tunneling problem 
through the impurity Green's function Ga while the po- 
sition dependence of the conductance is controlled by the 
electronic structure of the metal. 

(3) The spatial decay of the observed Fano resonance 
in the recent experiments^^'^^ is consistent with the con- 
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elusion that tunneling into the bulk conduetion and hy- 
bridized sp impurity states gives rise to most of the signal. 
The absence of any observable resonance at distances 
larger than ~ 10 A suggests that the contribution from 
the surface state on Au(lll) and Ag(lll) to the resonant 
tunneling is not important in these experiments. How- 
ever, the surface states are important in special cases, as 
indicated by the recent corral experiments^^ in which the 
contribution of the surface states is enhanced by scatter- 
ing from the walls of the corral. 

(4) At large Z^, tunneling into conduction states with 
having the smallest parallel component corresponding 

to energy oj = + efe^ is strongly favored. This leads 
to the disappearance of the current oscillation vs. the 
lateral tip position due to tunneling into the bulk states 
which should otherwise be observed with period of about 
1 — 2 A (corresponding to the bulk kp) for typical exper- 
imental conditions. Therefore no oscillations in the line- 
shape should be observed on this length scale for typical 
tip-surface separation. The occurrence of an antireso- 
nance with tip position at certain neighboring sites pre- 
dicted by Schiller and Hershfield^^ has its origin in these 
oscillations. It is a result of a simplified model for the 
surface electronic structure and we believe is unphysical. 
The small current oscillations predicted by Kawasaka et 
al.^^ assume that the surface states are all-important in 
the spatial dependence of the resonance which seems to 
contradict the experimental results. We believe the sur- 
face state should be important at larger distances since 
on the (111) noble metal surfaces kp ~ 0.15 — 0.2 A~^ 
and the corresponding period of oscillations is about 20 A 
(as observed experimentally as Friedel oscillations). We 
expect changes in the resonance lineshapes with this spa- 
tial period if the contribution from the surface state if the 
signature of the resonance is detectable at such distances. 

(5) From the lineshapes observed in Co/Au{lll) and 
Ce/Ag{lll), we conclude that the direct tunneling into 
the discrete (d or /) state is quite weak - stronger in 
Co/Au{lll). This confirms that the STM is mostly a 
probe of the delocalized sp states and couples only weakly 
to the tightly bound dor f orbitals at typical tip-surface 
separations. Therefore the dominant process giving rise 
to the resonance lineshape is the tip-to-metal tunneling 
and interference between conduction electrons scattering 
from the local moment. 

(6) The temperature dependence in differential con- 
ductance does not reflect only the temperature depen- 
dence of the Kondo resonance, but includes also the ef- 
fect of Fermi surface broadening (mostly of the tip). The 
two contributions are of the same order of magnitude and 
qualitatively indistinguishable. Therefore, the tempera- 
ture dependence in the differential conductance cannot 
be used directly to make conclusions about the tempera- 
ture dependence of the resonance without controlling the 
tip Fermi surface broadening or without deconvolution. 

(7) At small tip-surface separations, nonequilibrium ef- 
fects as well as the additional tip-adsorbate hybridiza- 
tion may play an important role - especially in Kondo 



systems. The main effect of the finite bias voltage in 
this case is to broaden the Kondo resonance and pro- 
duce nonequilibrium electron population on the adsor- 
bate. The observed Fano resonance in differential con- 
ductance also broadens and its maximum is suppressed. 
The effect of the tunneling current on the Kondo reso- 
nance should thus leave a characteristic dependence of 
the lineshape on Zt- 
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Appendix A. ADSORBATE GREEN'S FUNCTION 

An important quantity in the theory of tunneling 

current through adsorbate resonances is the adsorbate 
Green's function Ga- Using the equation of motion 
method, we find the expression for Ga defined as the 
Fourier transform of 

Ga{t,t') = -i{TcCa{t)ci{t')). (Al) 

We do this for the case of arbitrarily strong coupling 
between the tip and the adsorbate with the intent to de- 
scribe the nonequilibrium effects at finite bias. However, 
in this paper we consider the effect of the direct tip-metal 
interaction on Ga to be weak and neglect it. Extension 
to the full description will be considered in future work. 
We believe the approximations adopted here capture the 
most important nonequilibrium effects. 

We discuss both the noninteracting {U = 0) and inter- 
acting {U = oo) model. Since the solution in both lim- 
its for the adsorbate-metal interaction is well known, we 
limit our discussion to the issues specific to the addition 
of the biased tip and refer reader to standard texts for 
the details. The {U = oo) model is solved using the slave 
boson technique and NCA. In this approach a new pseud- 
ofermion is introduced by the transformation Ca — > Cab^ 
in the Hamiltonian (2.1), where is the creation oper- 
ator for the slave boson. This eliminates the interaction 
term U from the Hamiltonian as discussed by Coleman.®^ 

The time ordering operator Tc orders the time ac- 
cording to their position on contour in the complex time 
plane'^''. It is important to note that the equations must 
be first solved in the complex time domain and then an- 
alytically continued to the real axis as was pointed out 
in the previous appendix. The analytic continuation is 
performed before the Fourier transform, so we must be 
careful about how we deal with the Fourier transformed 
equations. Relevant details are in Appendix C. Here we 
discuss the equations of motion satisfied by the Fourier 
transforms of the time ordered Green's functions and the 
analytic continuation is performed at the end according 



24 



to the rules in Appendix C. All Green's functions and 
self-energies in the following expressions are function of 
frequency uo and, therefore, we omit their argument to 
simplify the notation. 

Gkk' 



Gpp' Gp 

► = ► + 



tpk hp' 




Vak Vk'a 
>c- *% I + 



'■ap '^p'a 

FIG. 13. (a) Diagrammatic expansion for the metal (G^fc/) 
and tip (Gpp/) Green's functions in mutual interaction but 
without the adsorbate. (b) diagrammatic expansion for the 
adsorbate Green's function Ga of the noninteracting Anderson 
model using the solutions G^k' and Gp^i . 

The Green's function for the impurity state Ga can be 
written in a standard way 

G„ = (w - eo - (A2) 

using the self energy T,a{Rt, Zq;uj). The solution for Ga is 
thus reduced to finding We first treat a closed shell or 
nonmagnetic open shell (Vka ^ U) adsorbate for which 
electron correlations can be neglected. We begin by con- 
sidering the tip-substrate system without the adsorbate. 
We define Gkk' and Gpp' in analogy with Ga (Al) as 
the Green's functions of the metal and tip states, re- 
spectively, in the absence of the adsorbate. These are 
not identical with the Green's functions Gkk' and Gppi 
for the full system introduced in Appendix B. The bare 
metal-tip system is described by the Hamiltonian of sec- 
tion II with f-o = U = tap = Vak = 0. Using the equations 
of motion, we can write 

{(J - Ck) Gkk' = Skk' + ^2 ^kk"Gk"k' (A3) 

k" 

and 

{u — Cp) Gpp' = Spp' + ^ ^ YippiiGpiipi (A-4) 
p" 



where the self-energies arc Y^kw = X]p tkpGptpk' and 

^PP' = Hk^pk^k^kp' and = {uj - Ck + irik)~^ and 
Gp = (u — ep + ir]p)~^ are the Green's functions for 
the clean metal and tip, respectively, without their mu- 
tual interaction. The coupled equations are shown di- 
agrammatically in Fig. 13. The solutions for Gkk' and 
Gppi can be formally written as the inverse of Dkk' = 
Skk'{'^-ek)-^kk' andDpp' = 6pp' {lu - ep) -'Epp' . We also 
define the adsorbate-metal and adsorbate-tip hybridiza- 
tion matrices modified by the tip-substrate interaction as 

Vka = Vka + tkpGptpa and ipa = tpa + X^fe tpkG'j^Vka- 

With these definitions and with Gkk' and Gpp' obtained 
through (A3) and (A4), the solution for the noninteract- 
ing Eq - shown diagrammatically in Fig. 13 - is formally 
given by 

= ^ VakGkk'Vk'a + ^ tapGpp'tp'a- (A5) 
kk' pp' 

The evaluation of the self energy Ea is rather complicated 
in the general case of strong tip-to-substratc coupling. 
We proceed with formulation of the general nonequi- 
librium theory for the tunneling current using this self 
energy (section III A) and then we discuss two limit- 
ing cases: (a) the equilibrium limit |tfep|,|iap| ^ \Vka\ 
(section IIIC) in which case the second term in (A5) 

is neglected and Vka, Gkk' replaced by Vka, G^Skk'', 
and (b) the nonequilibrium case under the assumption 
\tkp\ \tap\ ^ \Vka\, in which case we keep both terms 
in (A5) and replace 14a, tap, Gkk', and Gpp' by Vka, tap, 
GkSkk', and G^Spp'. Section IV B 3 deals with tunnel- 
ing through a Kondo impurity in this limit. The case 
(b) includes the effect of the increased hybridization of 
the discrete state due to the tip presence and the onset 
of nonequilibrium population on the adsorbate at finite 
bias. 

In order to study these corrections in the limit (b), 
{\tpk\ ^ \tpa\ ~ l^ftal)) we replace the Green's functions 
Gkk' and Gpp' by the noninteracting ones, i.e. Gkk' = 

Skk'Gl and Gpp' = Spp'Gp and the modified Vka and tpa 
by Vka and tpa- The self energy Eq then simplifies to 

S°=El^«fcl'G°fe + ElVl'G°. (A6) 
fc p 

The largest source of error in writing the approximate 
self energy is the neglect of the possibly significant in- 
terference efi'ects at larger tkp as a result of the phase 
difference between ipa and tpa- It is always reasonable 
to replace Vm by Vka , as long as the adsorbate is on the 
surface rather than on the STM tip. These general case 
will be the topic of a future study. If the tip distance 
from the adsorbate is much larger than the adsorbate- 
metal separation, so that tap <S Vak, the self-energy is 
well described by the first term only. In such a case, the 
STM docs not strongly modify the studied system. It is 
then reasonable to characterize the system without the 
presence of the STM tip and then consider the tunneling. 
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Finally, wc discuss the Green's function Ga in the limit 
(b) for a Kondo impurity which is likely to show stronger 
dependence on the bias and tip interaction. We find the 
nonequilibrium Green's function Ga under the same as- 
sumption that lead to E° for the noninteracting Ander- 
son Hamiltonian. We solve the interacting system in the 
limit oi U = oo using the NCA approximation, shown 
diagrammatically in Fig. 14. The self energy is not a 
simple sum of the two contributions from the metal and 
tip as it was in the noninteracting system, because the oc- 
cupation of the resonance is limited to one electron and 
the hybridization is now correlated formally through 
the slave boson Green's fimction B{u)). The two cou- 
pled equations in Fig. 14 are solved selfconsistently. The 
expressions for the NCA self energy is obtained in a stan- 
dard with the help of the diagrams in Fig. 14.^^ 



Ga 



onance and modifies the spectroscopic properties of the 
system. We therefore proceed by deriving the most gen- 
eral expression valid for arbitrary coupling strength tkp 
and tap and discuss an approximation (b) that allows us 
to take into account the most important nonequilibrium 
effects as described in the previous appendix. For this 
purpose we introduce Green's function Gpa{io), Gpk{co), 
Gkk'i'jj), Gpp'{uj), and Gka{i^) as the Fourier transform 
of 



Gij{t,t') = -i{TcCi{t)cUt')). 



(Bl) 



We now turn to the equations of motion for the Green's 
functions relevant for the tunneling current. The follow- 
ing expressions are valid for arbitrary interaction U 
and the nature of the intra-adsorbate interactions are 
contained fully in the solution for Ga discussed in the 
previous appendix. The first term in the current (3.3) 
contains the tip-adsorbate propagator which satisfies 




ap f 

-AAA*! 



t 

i pa 



(w — Cp) Gpa — tpaGa 



k 



tpkGka . 



(B2) 



It is expressed in term of Ga already solved within a given 
approximation in the previous appendix through (A2) 
and in terms of the metal-adsorbate Green's function 



(w — Ck) Gka = VkaGa + ^ tkpG 



pa' 



(B3) 



The last two equations are coupled and need to be solved 
self-consistently. We do this by substituting Eq. (B3) 
for G ka in (B2), and vice versa. The solutions are then 
expressed in terms of Ga, Gppi 
previous appendix as 



and tpa discussed in the 



pi 



ip'aGa' 



(B4) 



FIG. 14. The diagrammatic expansion for the adsorbate 
Green's function Ga of the interacting Anderson Hamiltonian 
{U = oo) within the NGA. 



Appendix B. EQUATIONS OF MOTION FOR Gpk 
AND Gpa 



In this appendix, we find the solution for 
(Y^kp^kpGpk + J2ap*apGpa) entering the expression for 
tunneling current (3.3) for the general case of arbitrary 
tip-system coupling. Ultimately, the interesting regime in 
connection with typical STM experiments is one in which 
tkpitap < Vak- Howcvci, wc waut to bc able, in principle, 
to study the system when the coupling of the STM tip 
to the system and the tunneling current are strong. This 
creates nonequilibrium occupation on the adsorbate res- 



We will also need the solution for Gak ■ The tip-induced 
correction to Vka contributes to the phase of Vka , as well 
as its magnitude, and could thus affect the lineshape sig- 
nificantly in the strong coupling limit. But it should be 
particularly weak when tkp, tpa <C Vka and it will be safe 
to ignore it. We write 



Gak — Ga Vak'Gk'k- 



(B5) 



The second term in (B3) is negligible when the tip- 
adsorbate separation is much larger than the adsorbate- 
metal separation. Neglecting this term is equivalent to 
replacing Gpp' — > G'^Spp' in (B4) and Gkk' G\5kk', 
Vka Vka in (B5). The tip-adsorbate Green's function 
Gpa is then expressed entirely in terms of Ga and the 
unperturbed conduction electron Green's functions. 

The second term in (3.3) contains the tip- metal prop- 
agator Gpk, which satisfies 
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(w - Cp) Gpk = '^tpaGak + '^tpk'Gk'k- (B6) 

a k' 

It is expressed in terms of Gak (B5) discussed in the 
previous paragraph and in terms of Gk'k, the Green's 
function for the substrate conduction electrons 

{UJ - Ck') Gk'k = hk' + ^ Vk'aGak + ^ h'pGpk- (B7) 
a p 

We see that Gkk' couples to Gpk, (B6), and also to 
Gak, (B5), already solved in terms of Go and Gkk'- The 
last two equations can be solved self-consistently to give 

Gk'k = Gk'k + ^ Gk'kiVakiGaVk2aGk2k- (B8) 
felfeg 

and 

Gpk = Gpp' {tp'k + ip'aGaVaki) G^. (B9) 
p' ak' 

For the purpose of analytic continuation, it is important 

to keep track of the order in which the Green's fimc- 
tions appear in the product in the above equations. The 
"lesser" Green's functions are then obtained according 
to rules stated in Appendix C. The equilibrium limit of 
the theory is achieved by neglecting the last term in (B7) 
along with the equivalent approximations for Ga and Gak 
discussed above. This removes the self-consistency re- 
quirement and neglects the effect of the tip on the sub- 
strate conduction electrons, but not on the tunneling cur- 
rent. The solution for Gkk' is then identical to that of 
the system without the tip. 

Appendix C. RULES FOR DEALING WITH 
NONEQUILIBRIUM GREEN'S FUNCTIONS IN 
FREQUENCY SPACE 



/OO 
rfri...rfr2B«(^)(t,ri)...Z«(^)(r2,0 
-OO 

(C2) 

which consists of only one term and 

/OO 
dTi...dTs...dT2 (C3) 
-OO 

[... + S^(t,Ti)...C^ {t,rs)...Z^{r2,t') + ...] 

where each of the N terms in the integral has exactly one 
function of the type / ^ , all functions to the left (right) 
of it arc retarded (advanced) , and each of the terms has 
the / ^ in a different position. 

We are dealing here with the case of time independent 
perturbations (steady state current). The double time 
propagators then only depend on the time difference {t — 
t'). The equations (C2), (C3) can be Fourier transformed 
and we can write the time ordered Green's function A{lo) 
as a simple product 

A{oj) = B{lo)...Z{oj). (C4) 

The rules for writing the expression for the "lesser" and 
retarded function in the frequency space are then directly 

carried out by Fourier transforming the equations (C2) 
and (C3). Leaving out the frequency arguments, we write 

and 

A^ = ... + B^...C^ ...Z^ + ... (06) 

We note that it is important to keep track of the order 
in which the function B...Z appear in the time integral 
when the Fourier product is formed. 



In this appendix we review the process of analytic con- 
tinuation of the complex time contour expression to inte- 
grals on the real time axis. We follow Langreth's general- 
ization^^ of Kadanoff-Baym's method"^^ described in de- 
tail by Haug and Jauho^^. Four Green's functions appear 
in the nonequilibrium theory "lesser" G*^, "greater" G^, 
retarded G^, and advanced g?^ 46 48,26 _ frequently 
need to find the retarded (advanced) and "lesser" Green's 
function corresponding to Green's fimction A time or- 
dered in the complex time plane expressed as a product 
of N time ordered functions B...Z , i.e. 

A{t,t')= I dTi...dT2B{t,Ti)..Z{T2,t') (CI) 

Jc 

where all functions are assumed fermion-like. The desired 
expressions analytically continued onto the real time axis 

47 
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